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, Abstract 

| |The Cherenkov Telescope Array (CTA) is a project for a next-generation observatory for very high energy (GeV-TeV) ground-based 
'gamma-ray astronomy, currently in its design phase, and foreseen to be operative a few years from now. Several tens of telescopes 
of 2-3 different sizes, distributed over a large area, will allow for a sensitivity about a factor 10 better than current instruments 
such as H.E.S.S, MAGIC and VERITAS, an energy coverage from a few tens of GeV to several tens of TeV, and a field of view of 
up to 10 deg. In the following study, we investigate the prospects for CTA to study several science questions that can profoundly 
influence our current knowledge of fundamental physics. Based on conservative assumptions for the performance of the different 
CTA telescope configurations currently under discussion, we employ a Monte Carlo based approach to evaluate the prospects for 
detection and characterisation of new physics with the array. 

First, we discuss CTA prospects for cold dark matter searches, following different observational strategies: in dwarf satellite galaxies 
• • of the Milky Way, which are virtually void of astrophysical background and have a relatively well known dark matter density; in 
. £h \ < ^ ie re gi° n close to the Galactic Centre, where the dark matter density is expected to be large while the astrophysical background 
-due to the Galactic Centre can be excluded; and in clusters of galaxies, where the intrinsic flux may be boosted significantly by the 
$_i large number of halo substructures. The possible search for spatial signatures, facilitated by the larger field of view of CTA, is also 
. F^. discussed. Next we consider searches for axion-like particles which, besides being possible candidates for dark matter may also 
explain the unexpectedly low absorption by extragalactic background light of gamma-rays from very distant blazars. We establish 
the axion mass range CTA could probe through observation of long-lasting flares in distant sources. Simulated light-curves of 
flaring sources are also used to determine the sensitivity to violations of Lorentz Invariance by detection of the possible delay 
between the arrival times of photons at different energies. Finally, we mention searches for other exotic physics with CTA. 

Keywords: CTA, Dark Matter, Dwarf satellite galaxies, Galactic centre, Galactic halo, Galaxy clusters, Axion-like Particles, 
Lorentz Invariance Violations, Neutrino, Magnetic monopoles, Gravitational Waves 
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Introduction 

The Cherenkov Telescope Array (CTA) [1] will be an ad- 
vanced facility for ground-based gamma-ray astronomy in the 
GeV-TeV regime. Compared to the current generation of Imag- 
ing Atmospheric Cherenkov Telescopes (IACT) e.g. H.E.S.S., 
MAGIC and VERITAS 2 , CTA will feature substantial improve- 
ments. It will cover over 3 decades in energy, from a few tens of 
GeV up to several tens of Te V. At both ends of this range, inter- 
esting new physics is expected, and in addition, a larger energy 
coverage will provide a bigger lever arm for spectral studies. 
Above 1 TeV, the field of view (FOV) will be up to 10 deg 
i.e. over a factor of 2 larger than that of current instruments. 
CTA is currently planned to have a Southern hemisphere site 
and Northern hemisphere site. This fact together with the large 
FOV of the telescopes in both installations will likely enable 
CTA to provide the first extended gamma-ray maps of the sky 
in the TeV region. The improved energy and angular resolution 
will enable more precise spectral and morphological observa- 
tion. This will be achieved by deploying several tens of tele- 
scopes of 2-3 different sizes over an area of several square km. 
CTA will be operated as an open observatory, with improved 
data dissemination among the world-wide scientific community 
and a substantial fraction of the total observation time devoted 
to guest proposals. 

The search for new physics beyond the Standard Model 
(SM) of particle physics is among the key science drivers of 
CTA along with the understanding of the origin of high-energy 
gamma-rays and of the physics of cosmic ray acceleration in 
galactic and extragalactic objects. Several such fundamental 
physics issues are examined here — the nature of cold dark 
matter, the possible existence of axion-like particles, and ex- 
pected violation of Lorentz Invariance by quantum gravity ef- 
fects. Search strategies for cosmic tau neutrinos, magnetic 
monopoles and follow-up observations of gravitational waves, 
are also discussed. 

The CTA array performance files and analysis algorithms 
are extensively described in Bernloher et al. [2]. Eleven 
array configurations (A...K) were tested for the Southern 
hemisphere and two (NA, NB) for the Northern hemisphere 
[2, Table 2]. The simulations were made at an altitude of 2000 
m and at 70 deg elevation. Arrays E and / are considered 
balanced layouts in terms of performance across the energy 
range. Arrays A, B, F and G are more focused to low-energies, 
and arrays C, D and H to high energies. NB is a higher 
energy alternative to NA. Their point-source sensitivity is 
compared in [2, Fig. 7]. The arrays comprise different number 
of telescopes of three different sizes: the Large Size Telescope 
(LST, 23 m diameter), the Medium Size Telescope (MST, 
12 m diameter) and Small Size Telescope (SST, 6 m diameter) 
[2, Table 1]. One of the goals of this study was to compare 
different array configurations for the specific scientific case. 
While in some cases all CTA configurations are compared 
against each other, in others only benchmarks array B, C and 
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E are considered, as representative arrays that maximize the 
performance at low-energy, high-energy and in the full-range, 
respectively. Except for galaxy cluster studies and Galactic 
halo studies, where extended or diffuse MC simulations are 
used, in all other cases point-like MC simulations are used. 
This is the first time that realistic estimates of the prospects 
of detection for CTA are presented for such searches. An 
optimised event selection procedure and a dedicated analysis 
ought to improve on our conservative expectations. Previous 
studies often relied on too optimistic sensitivities, especially 
at low energies (< 100 GeV); publicly available effective 
areas for a subset of configurations [1, 2] are now accurate 
and can be used to infer CTA sensitivities for point-like sources. 

This contribution is structured as follows: 

- In Section 1, we explore different possible scenarios for 
detection of cold dark matter particle signatures in obser- 
vations of: dwarf satellite galaxies of the Milky Way (Sec- 
tion 1.1), clusters of galaxies (Section 1 .2) and the Galactic 
halo (Section 1.3). We also study anisotropies in the dif- 
fuse gamma-ray background as a signature of dark matter 
(Section 1.4). 

- In Section 2, we discuss the scientific case for axion-like 
particles, and make predictions for detection from obser- 
vation of blazars at different distances and with different 
flare durations. 

- In Section 3, we compare the capacity of all planned CTA 
arrays to constrain high energy violations of Lorentz In- 
variance, relative to current limits. 

- In Section 4 we discuss qualitatively three more cases: 
the observation of air showers from T-leptons emerging 
from the Earth's crust (Section 4.1), the capability to iden- 
tify magnetic monopoles as bright emitters of Cherenkov 
light in the atmosphere (Section 4.2) and some considera- 
tion about multi-wavelength gravitational wave campaigns 
(Section 4.3). 

Given the wide variety of physics issues considered in this 
contribution, an introduction to the individual physics case is 
presented in each section for easier readability. The reader can 
find an overall summary and closing remarks in Section 5. 

1. Cold Dark Matter Particle searches 

A major open question for modern physics is the nature of 
the dark matter (DM). There is a large body of evidence for the 
presence of an unknown form of gravitational mass, at scales 
from kiloparsecs to megaparsecs, that cannot be accounted for 
by SM particles. The observation by the WMAP satellite [3] 
of the acoustic oscillations imprinted in the cosmic microwave 
background quantifies the DM component as contributing about 
25% of the total energy budget of the Universe. Being dominant 
with respect to the baryonic component, which accounts for 
only about 4% of the total energy density, DM shaped the for- 
mation of cosmic structures. By comparing the galaxy distribu- 
tions in large redshift galaxy surveys [4], and through /V-body 
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simulations of structure formation [5, 6, 7], it is inferred that 
the particles constituting the cosmological DM had to be mov- 
ing non-relativistically at decoupling from thermal equilibrium 
in the early universe ( 'freeze-out' ), in order to reproduce the ob- 
served large-scale structure in the Universe and hence the term 
"cold DM" (CDM). This observational evidence has led to the 
establishment of a concordance cosmological model, dubbed 
ACDM [8, 9, 10], although this paradigm is troubled by some 
experimental controversies [11, 12, 13, 14, 15, 16]. 

One of the most popular scenarios for CDM is that of weakly 
interacting massive particles (WIMPs), which includes a large 
class of non-baryonic candidates with mass typically between 
a few tens of GeV and few TeV and an annihilation cross- 
section set by weak interactions [see, e.g., Refs. 17, 18]. Nat- 
ural WIMP candidates are found in proposed extensions of the 
SM, e.g. in Super- Symmetry (SUSY) [19, 20], but also Lit- 
tle Higgs [21], Universal Extra Dimensions [22], and Tech- 
nicolor models [23, 24], among others. Their present veloci- 
ties are set by the gravitational potential in the Galactic halo at 
about a thousandth of the speed of light. WIMPs which were 
in thermal equilibrium in the early Universe would have a relic 
abundance varying inversely as their velocity-weighted annihi- 
lation cross-section (for pure s-wave annihilation): ficDM^ 2 = 
3 x 10 _27 cm 3 s -1 / (cr ann v) [19]. Hence for a weak-scale cross- 
section ((Tannv) = 3 x 10 _26 cm 3 s _1 , they naturally have the 
required relic density Ocdm^ 2 = 0.113 ± 0.004, where h - 
0.704 ± 0.014 is the Hubble parameter in units of 100 km s -1 
Mpc 1 [3]. The ability of WIMPs to naturally yield the DM 
density from readily computed thermal processes in the early 
Universe without much fine tuning is sometimes termed the 
"WIMP miracle". 

In some SUSY theories, a symmetry called '/^-parity' 
prevents a too rapid proton-decay, and as a side-effect, also 
guarantees the stability of the lightest SUSY particle (LSP), 
which is thus a prime candidate for a WIMP. WIMPs can 
annihilate to SM particles, and have hadron or leptons in 
the final products of annihilation. Thus from cosmic DM 
annihilations, one can expect emission of neutrinos, charged 
cosmic rays, multi-frequency electromagnetic radiation from 
charged products, and prompt gamma-rays [25]. The detection 
of these final state particles can help to identify DM — this is 
termed "indirect DM detection". Gamma-rays are not deflected 
by cosmic magnetic fields, and thus trace back to their origin. 
Therefore, observation of a gamma-ray signal from cosmic 
targets where DM is expected could prove conclusive about its 
nature . 

In the context of gamma-ray astronomy, the differential flux 
of gamma-rays from within a solid angle AO around a given 
astronomical target where DM is expected, can be written as: 

«. / Astrophysics 

Particle Physics 

(l.i) 

where (<x ann v) is the annihilation cross-section (times the rela- 



tive velocity of the two WIMPs), BR; dA^/dE r = dN 7 /dE y 
is the photon flux per annihilation summed over all the possible 
annihilation channels i with branching ratios BR,, and m x is the 
mass of the DM particle. The 'astrophysical factor' J is the in- 
tegral over the line of sight (los) of the squared DM density and 
over the integration solid angle AO: 

7= f ffif Asp 2 (s,Q.). (1.2) 

JACl Jlos 

The remaining term B F in Eq. (1.1) is the so-called 'boost fac- 
tor' which is a measure of our ignorance of intrinsic flux con- 
tributions that are not accounted for directly in the formula. 

There are various known mechanisms for boosting the intrin- 
sic flux, among which we mention the inclusion of subhalos, 
and the existence of a 'Sommerfeld enhancement' of the cross- 
section at low velocity regimes in models where the DM parti- 
cles interact via a new long-range force. All numerical A^-body 
simulations of galactic halos have shown the presence of sub- 
halos populating the host halo [see, e.g., Refs. 5, 26]. Such 
density enhancements, if not spatially resolved, can contribute 
substantially to the expected gamma-ray flux from a given ob- 
ject. This effect is strongly dependent on the target: in dwarf 
spheroidal galaxies (dSphs) for example the boost factor is only 
of 0(1) [27, 28], whereas in galaxy clusters the boost can be 
spectacular, by up to a factor of several hundreds [29, 30, 31]. 
On the other hand, the Sommerfeld enhancement effect can 
significantly boost the DM annihilation cross-section [32, 33]. 
This non-relativistic effect arises when two DM particles inter- 
act in a long-range attractive potential, and results in a boost 
in gamma-ray flux which increases with decreasing relative ve- 
locity down to a saturation point which depends on the DM and 
mediator particle mass. This effect can enhance the annihilation 
cross-section by a few orders of magnitude [27, 28]. 

The current generation of IACTs is actively searching for 
WIMP annihilation signals. dSphs are promising targets for 
DM annihilation detection being among the most DM domi- 
nated objects known and free from astrophysical background. 
Constraints on WIMP annihilation signals from dSphs have 
been reported towards Sagittarius, Canis Major, Sculptor and 
Carina by H.E.S.S. [34, 35, 28], towards Draco, Willman 1 and 
Segue 1 by MAGIC [36, 37, 38], towards Draco, Ursa Minor, 
Bootes 1, Willman 1 and Segue 1 by VERITAS [39, 40], 
and again towards Draco and Ursa Minor by Whipple [41]. 
Nevertheless, the present instruments do not have the required 
sensitivity to reach the "thermal" value of the annihilation 
cross-section (<x ann v) = 3x 10 _26 cm 3 s _1 . A search for a WIMP 
annihilation signal from the halo at angular distances between 
0.3° and 1.0° from the Galactic Centre has also recently been 
performed using 112h of H.E.S.S. data [42]. For WIMP 
masses well above the H.E.S.S. energy threshold of 100 GeV, 
this analysis provides the currently most constraining limits 
on (o- ann v) at the level of a fewxl0~ 25 cmV 1 . H.E.S.S., 
MAGIC and VERITAS have also observed some galaxy 
clusters, reporting detection of individual galaxies in the 
cluster, but only upper limits on any CR and DM associated 
emission [43, 44, 45, 46, 47, 48]. Even though IACT limits 
are weaker than those obtained from the Fermi-LAT satellite 
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measurements in the GeV mass range [49, 50, 51, 52], they 
complement the latter in the TeV mass range. Gamma-ray line 
signatures can also be expected in the annihilation or decay of 
DM particles in space, e.g. into yy or Z°y. Such a signal would 
be readily distinguishable from astrophysical gamma-ray 
sources which typically produce continuous spectra [53]. A 
measurement carried out by H.E.S.S. [54] using over 100 h of 
Galactic Centre observations and over 1000 h of extragalactic 
observations complements recent results obtained by Fermi- 
LAT [55], and together cover about 3 orders of magnitude in 
energy, from 10 GeV to 10 TeV. 

In this contribution, we focus on the prospects for DM 
searches with CTA, which are expected to improve on the cur- 
rent generation of IACTs on the following basis: 

- the energy range will be extended, from a few tens of GeV 
to several tens of TeV. At low energies, this will allow 
overlap with the Fermi-LAT instrument, and will provide 
sensitivity to WIMPs with low masses. For WIMPs with 
mass larger than about 100 GeV, CTA will have higher 
sensitivity as our studies indicate [56]. 

- the improved sensitivity in the entire energy range, com- 
pared to current instruments, will obviously improve the 
probability of detection, or even identification of DM, 
through the observation of spectral features, 

- the increased FOV (about 10 deg versus 2-5 deg) with 
a much more homogeneous sensitivity, as well as the im- 
proved angular resolution, will allow for much more ef- 
ficient searches for extended sources like galaxy clusters 
(Section 1 .2) and spatial anisotropics (Section 1 .4), 

- finally, the improved energy resolution will allow much 
better sensitivity to the possible spectral feature in the DM- 
generated photon spectrum. While astrophysical sources 
show typically power-law spectra with steepening at high 
energies, DM spectra are universal and generically ex- 
hibit a rapid cut-off at the DM mass. For specific mod- 
els, "smoking gun" spectral features can appear [53]. The 
observation of a few identical such spectra from differ- 
ent sources will allow both precision determination of the 
mass of the WIMP and its annihilation cross-section. 

For the following studies, in order to have a detection, we re- 
quire a) the number of excess events over the background larger 
than 10 in the signal region, b) the ratio between the number of 
excess events and the background events larger than 3%, and c) 
the significance of the detection computed following Eq. (17) 
of Li&Ma [57], to be larger than 5. If not explicitly mentioned, 
we used a number of background-control regions set to 5 (a 
= 0.2 in the Li&Ma notation), which is a conservative choice, 
given the fact that the large FOV of CTA may allow for a < 0.2. 
In case of non detection within a certain observation time, we 
calculate integral upper limits following the methods described 
in Rolke et al. [58] (bounded profile likelihood ratio statistic 
with Gaussian background, and with a confidence level of 95% 
C.L) in all cases expect the Galactic halo case, where we use 
the method of Feldman and Cousins [59]. 
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Ursa Minor 


+44.8 


66 


2.2 x 10 18 


NFW 


[66] 


Draco 


+34.7 


87 


7.1 x 10 17 


NFW 


[66] 


Sculptor 


-83.2 


79 


8.9 x 10 17 


NFW 


[66] 








2.7 x 10 17 


ISO 


[68] 


Carina 


-22,2 


101 


2.8 x 10 17 


NFW 


[66] 


Segue 1 


+ 16.1 


23 


1.7 x 10 19 


Einasto 


[38] 


Willman 1 


+51.1 


38 


8.4 x 10 18 


NFW 


[39] 


Coma Berenices 


+23.6 


44 


3.9 x 10 18 


NFW 


[69] 



Table 1.1: Astrophysical factors for a selection of the most promising classical 
and ultra-faint dSphs. Dec. is the target declination, D the distance and J is 
defined as in Eq. 1.2. 



We study the effect of various annihilation spectra, assuming 
in turn 100% BR into a specific channel (bb, t + t~ or // + yU )■ 
The spectral shapes are obtained from different parameterisa- 
tion from the literature [60, 61, 62]. For the bb channel, which 
is used for comparison of different targets (see Fig. 5.1), this 
difference accounts for few percents (depending on the DM 
mass), which is substantially smaller than the uncertainties in, 
e.g., the astrophysical factor, and do not significantly alters the 
conclusions. 

1.1. Observations of dwarf satellite galaxies 

In the ACDM paradigm, galaxies such as ours are the result 
of a complex merger history and are expected to have extended 
halos of DM in accordance with observations. dSphs are satel- 
lites orbiting the Milky Way under its gravitational influence 
and are considered as privileged targets for DM searches for 
the following reasons: 

- the study of stellar dynamics shows that dSphs are among 
the most DM-dominated systems in the Universe, with 
mass-to-light ratio up to a few hundreds. In particular, the 
otherwise very uncertain astrophysical factor (Eq. 1 .2) can 
be constrained by dynamical arguments [63], 

- many of the dSphs lie within ~ 100 kpc of the Earth, 

- they have favourable low gamma-ray backgrounds due to 
the lack of recent star formation history and little or no gas 
to serve as target material for cosmic -rays [64]. 

The family of dSphs is divided into "classical" dSphs, which 
are well-established sources with relatively high surface bright- 
ness and hundreds of member stars identified [65, 66], and 
"ultra-faint" dSphs, which have mainly been discovered re- 
cently through photometric observations in the Sloan Digital 
Sky Survey (SDSS) [67] and have very low surface brightness 
and only a few tens or hundreds of member stars. Some of 
the ultra-faint dSphs are not well-established as such because 
of similarity of their properties with globular clusters, hence 
their nature is often under debate. However, they are of partic- 
ular interest due to their potentially very large, albeit uncertain, 
mass-to-light ratios. 

Table 1 . 1 shows the astrophysical factor J for few selected 
dSphs for comparison. For the classical dSphs, we selected 
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the two most promising Northern (Ursa Minor and Draco) and 
Southern (Sculptor and Carina) ones according to Charbonnier 
et al. [66, Table 2]. The statistical uncertainties on the as- 
trophisical factor are roughly one order of magntiude at 68% 
CL, slightly depending on the dSphs, and can be found in 
[66, Table 2]. For the ultra-faint dSphs, we include Segue 1, 
Willman 1 and Coma Berenices, which have the highest J- 
values (although their nature is still under debate, especially for 
Segue 1 [70, 71, 72, 73, 74, 75, 76], which makes the determi- 
nation of the astrophysical factor less accurate than for classical 
dSphs). We remark how the estimation of the astrophysical fac- 
tor is subject to uncertainties of either statistical origin or due to 
the different assumptions considered for its calculation. A sys- 
tematic study has been done for Sculptor, to estimate the effect 
of the profile shape and velocity anisotropy assumptions [68]. 
Another compilation of astrophysical factors for several dSphs 
can be found in [52]. 

For the subsequent discussion, we consider only three 
sources: Ursa Minor and Sculptor representative of classic 
dSphs and located in the Northern and Southern hemisphere 
respectively, and Segue 1 having the largest astrophysical fac- 
tor. 

Bounds on the annihilation cross-section 
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Figure 1.1: CTA sensitivities on the velocity-averaged annihilation cross- 
section as a function of the WIMP mass for 100 hours observation of Sculptor 
with the CTA array E (solid line), B (dashed line) and C (dashed-dotted line). 
Both the NFW (black line) and cored isothermal (ISO, red line) DM halo pro- 
files are shown, for an integration solid angle A£l = 1 X 10~ 5 sr. Annihilations 
are assumed to occur with 100% branching ratio into bb. 

Two kinds of radial profiles are generally used to model 
the DM distribution in dSphs: cusped and cored profiles [77]. 
While the former is motivated by numerical Af-body simu- 
lations, the latter seems to be more consistent with observa- 
tions [78], but the issue is still under debate [see, e.g., 79]. The 
standard cusped profile is the Navarro, Frenk & White form 



(NFW) [80], while more recently it has been shown that the 
Einasto profile [81] provides also a good fit to the subhalos in 
A-body simulations [5]. On the other hand, for systems of the 
size of dSphs, the possibility of centrally cored profiles has also 
been suggested [82, 83, 14]. In conclusion, observations of low 
surface brightness and dSphs [84, 85, 86] show that both cusped 
and cored profiles can accommodate their stellar dynamics. 

Fig. 1 . 1 shows the integral upper limits towards Sculptor, the 
best Southern candidate from Table 1.1, for which we consider 
both a cusped NFW [66] and a cored isothermal [28] profile. 
The sensitivity is calculated assuming that the DM particle an- 
nihilates purely in the bb channel, for arrays B, C and E. The 
observation time is set to 100 hours and the integration solid 
angle to AO = 10~ 5 sr. The best reached sensitivity is at the or- 
der of fewxl0~ 23 cm 3 s -1 for the NFW profile for both arrays E 
and B, while the isothermal profile is less constraining. Weaker 
constraints in the low mass range are obtained for the C array 
due to the lack of the large-size telescopes in the centre of their 
layout. The capability of CTA to discriminate between the two 
profiles is therefore restricted. 

The integration solid angle plays a central role in the 
estimation of the sensitivity and in the discrimination of the 
cusp or core profiles. The former point was addressed already 
[66, Fig. 7] where it was shown that small integration angles 
guarantee the strongest constraints. In the case of CTA, 
depending on the array layout (and the energy range), the 
angular resolution could be as low as 0.02 deg, corresponding 
to a minimum integration angle of about 10~ 6 sr, and thus 
our results can be considered conservative, with an expected 
improvement of up to a factor ~ 2. Concerning the second 
point, Walker et al. [87] showed that the more robust con- 
straints, regardless of whether the profile is cored or cusped, 
are reached for an integration angle r c = 2ri/2/A where 
ri/2 is the so-called half-light radius, and D is the distance to 
the dSph. For Sculptor, r c - 0.52, which is over 5 times the 
integration angle adopted here. In our calculation this would 
imply a weakening of the upper limits of a factor of a few. 

In Fig. 1.2 we show the integral upper limits for two classi- 
cal dSphs, namely Ursa Minor and Sculptor in the Northern 
and Southern hemispheres respectively, as well as the ultra- 
faint dSph Segue 1. In order to span the variety of DM par- 
ticle models, we study the effect of various annihilation spectra 
(computed using Ref. [61]), assuming in turn 100% BR into 
bb, t + t~ and pfyT channels for the array E and an observa- 
tion time T b s = 100 h. Assuming the annihilation to be purely 
into t + t~, the sensitivity reaches few xl0~ 25 cm 3 s _1 for 100 h 
observation time of Segue 1 . In comparing the different dSphs 
(assuming the reference annihilation channel bb) we see that 
even the most promising classical dSphs are less constraining 
than Segue 1 by over a factor of 10. However the uncertainties 
in the estimation of astrophysical factors for ultra-faint dSphs 
mean that this conclusion may not be reliable. Note that in the 
above calculations we did not assume any intrinsic flux boost 
factor, i.e. Bp = 1 in Eq. (1.1). 
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Figure 1.2: CTA sensitivities on the velocity-averaged annihilation cross- 
section versus the WIMP mass for 100 hours observation towards Sculp- 
tor, Ursa Minor and Segue 1, assuming 100% branching ratio into bb (for 
Segue 1 also into r + r~ and The calculations are done for array E and 

AH = 1 X 10~ 5 sr. 
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Figure 1.3: The minimum value of the astrophysical factor required for a 5 cr 
detection with T [, s = 100 h, versus WIMP mass. Two annihilation channels are 
considered for arrays B, C, and E: bb (upper curves) and t + t~ (lower curves). 
The estimated astrophysical factor for Segue 1 is shown for comparison. 



Bounds on Astrophysical factors and Boost factors 

Another approach to estimate the capabilities of CTA for DM 
detection in dSphs consists in the evaluation of the statistical 
significance of the DM signal as a function of the DM particle 
mass m x and the astrophysical factor, for different possible 



annihilation channels. Hereafter, we calculate the minimum 
astrophysical factor 7 m ; n required to reach a statistical signifi- 
cance of 5 cr assuming an effective observation time of 100 h, 
and the thermal cross-section 3 x 10~ 26 cm 3 s~'. This is shown 
in Fig. 1.3 for two annihilation channels: bb (upper curves) 
and t + t~ (lower curves), using analytical fits from Ref. [88]. 
Again, three proposed CTA configurations are studied: B, C, 
and E. In order to put these values into context, we note that 
the largest astrophysical factor / for known dSphs is that of 
Segue 1 at 1.7 x 10 19 GeV 2 cnr 5 [89]. From the figure we 
see that array B is the most constraining over the whole energy 
range. It is clear that for a detection, the astrophysical factor of 
the dSph needs to exceed 10 21 GeV 2 cm" 5 , which is only 1-2 
orders of magnitude smaller than that of the Galactic Centre 
(see Section 1.3). While we may expect a few such objects 
in the Milky Way halo [90], they ought to have already been 
detected and identified by Fermi-LAT. Although this has not 
happened, one can envisage DM subhalos with no associated 
dSph (or one not bright enough optically to be detected), and 
therefore such gamma-ray emitters may be hidden among the 
unidentified Fermi sources [91]. 

Another way to evaluate the prospects of DM detection is by 
means of the intrinsic flux boost factor term Bp in Eq. (1.1). 
The minimum Bp is computed as the ratio of the minimum as- 
trophysical factor / m j n which provides a 5 cr detection in 100 h 
of observation time with CTA, to the observational astrophysi- 
cal factor J from the DM modeling of the dSphs. Again, the 
thermal cross-section 3 x 10" 26 cm 3 s"' is assumed. Fig. 1.4 
shows the minimum Bp for a 1 TeV DM particle annihilating 
into t + t~. 7 b s is calculated for a NFW profile for all the cases 
except Segue 1, where an Einasto profile is considered. Con- 
sidering that the boost factor from subhalos in dSph is only of 
0(1), CTA observations of dSphs will be more sensitive to sce- 
narios where Sommerfeld enhancement is at play, which may 
instead boost the signal up to <9(1000). 
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Figure 1.4: Minimum boost factor required for a 5cr detection in 100 h by 
array B, for the dSphs in Table 1.1 and a 1 TeV WIMP annihilating into r + r~. 
The density profiles are taken to be NFW, except for Segue 1 where an Einasto 
profile has been assumed. The smallest boost required is Bp = 25 for Segue 1 . 
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1.2. Observations of Galaxy Clusters 

Within the standard ACDM scenario, galaxy clusters, with 
masses around 10 14 - 10 15 M Q , are the largest gravitation- 
ally bound objects and the most recent structures to form [92]. 
They are complex objects, relevant for both cosmological and 
astrophysical studies, and for what concerns DM searches 
[25, 29, 30, 93, 94, 95, 96, 97, 98, 99, 100]. DM, in fact, 
is supposed to be the dominant component of the cluster mass 
budget, accounting for up to 80% of its mass (the other compo- 
nents are the galaxies and the gas of the intra-cluster medium 
(ICM)). This is why clusters have been considered as targets for 
the indirect detection of DM, with the possibility of detecting 
the gamma-rays produced in the annihilation (or decay) of DM 
particles in the halo of the cluster. 

Af-body simulations of halo formation and evolution have 
also proven that, while the majority of early-formed, small 
structures merge together giving shape to more massive objects, 
some of the subhalos survive and are still present in the "host" 
halo of larger objects. Theoretical models foresee a huge num- 
ber of these substructures at all scales down to 10~ n - 10~ 3 M Q 
[101]. These subhalos have the effect of contributing to the to- 
tal gamma-ray emission from DM annihilations, and they may 
have important consequences for DM indirect detection. This 
is especially true for galaxy clusters, where the intrinsic flux 
"boost" from subhalos can be of order 100 - 1000, in particu- 
lar compared to the case of dSphs, explored previously, where 
the subhalos boost should contribute only marginally. Despite 
the fact that, due to their vicinity, dSphs are usually consid- 
ered as the best sources for DM indirect detection, thanks to the 
subhalos boost, some authors claim that galaxy clusters have 
prospects of DM detection better or at least as good as those of 
dSphs [29, 30,31]. 

On the other hand, in galaxy clusters, emission in the gamma- 
ray range is not only expected by DM annihilation. Clus- 
ters may host an Active Galaxy Nucleus (AGN, that appear 
as point-like sources at very high energies) and radio galax- 
ies. The case of the Perseus galaxy cluster, which has been 
observed by MAGIC during several campaigns in the last 
years, is emblematic: MAGIC detected both the central AGN 
NGC-1275 [102] and the off-centreed head-tail radio galaxy 
IC 310 [103]. Moreover gamma-rays are expected to be pro- 
duced also from the interaction of cosmic rays (CRs) with the 
ICM [95, 104, 105, 106, 107]. The physics of the accel- 
eration of CRs (electrons and protons) is not completely un- 
derstood, but plausible mechanisms can be shock acceleration 
during structure formation, or galactic winds driven by super- 
novae. CRs can also be injected into the ICM from radio galaxy 
jets/lobes. At the energies of interest here (above 10 GeV), 
CRs emit gamma-rays from the processes associated with the 
decay of the neutral and charged pions produced in the inter- 
action of the CRs with the ICM ambient protons [108, 109]. 
Most importantly, such a contribution is usually found to be 
larger than the one predicted from DM annihilation. It thus rep- 
resents an unavoidable source of background for DM searches 
in galaxy clusters. To date, the deep exposure performed with 
the MAGIC stereoscopic system of the Perseus cluster [47] 
placed the most stringent constraints from VHE gamma-rays 



observations regarding the maximum CRs-to-thermal pressure 
to {X CR ) < 1 - 2%. 

The purpose of this section is to estimate the CTA poten- 
tial to detect gamma-rays from DM annihilation in the halo of 
galaxy clusters. First, the CR-induced emission only will be 
considered. This component represents, by itself, an extremely 
interesting scientific case, at the same time being a background 
complicating the prospects of DM detection. Afterward, the 
ideal case of a cluster whose emission is dominated by DM an- 
nihilation only will be treated. Finally, the combination of the 
two components distributed co-spatially will be discussed. 

It should be noted here that gamma-ray emission from both 
DM annihilation and CRs is spatially extended, even though 
not always co-spatial. In particular, Sanchez-Conde et al. [29] 
proved that, for the case of DM, the contribution of subhalos is 
particularly relevant away from the halo centre, so that annihi- 
lations can still produce a significant amount of photons up to 
a distance of 1 - 2 degrees from the centre. This represents a 
problem for current Cherenkov Telescopes since their FOV is 
limited to 3 - 5 degrees. CTA will overcome this limitation, 
having a FOV of up to 10 deg (at least above 1 TeV) and an 
almost flat sensitivity up to several degrees from the centre. It 
is reasonable to expect, therefore, that CTA will allow a step- 
change in capability in this important area. 

In this study, we selected two benchmark galaxy clusters: 
Perseus and Fornax. Perseus has been chosen because it 
is considered that with the highest CR-induced photon yield 
but a low DM content, and Fornax for the opposite reason: 
it is considered the most promising galaxy cluster for DM 
searches [29, 30]. We recall that Perseus is located in the North- 
ern hemisphere, while Fornax is in the Southern hemisphere. 
To study the prospects for CTA we use two Monte Carlo sim- 
ulations of the instrument response functions and of the back- 
ground rates for extended sources, for the case of array B and 
array E, which we recall, are representatives of well-performing 
arrays at low energies (array B) and in the full energy range (ar- 
ray E). The MC simulations were developed explicitly for the 
analysis of extended sources so that all the relevant observables 
are computed throughout the entire FOV. 

Gamma-ray emission from cosmic-rays 

Gamma-ray emission due to the injection of CRs into the 
ICM of a galaxy cluster is proportional both to the density of 
the ICM and the density of CRs. For the present work, we refer 
to the hadronic CR model of Pinzke and collaborators [30, 95], 
based on detailed hydrodynamic, high-resolution simulations of 
the evolution of galaxy clusters, since in these works we found 
detailed morphological information, essential to compute the 
CTA response. The CR surface brightness rapidly decreases 
with the distance from the centre of the halo, so that, in most 
cases, the total emission is contained in 0.2 - 0.3 r 2 oo, where 
r2oo is the projected virial radius of the cluster, where the local 
density equals 200 times the critical density (see, e.g., Fig. 14 of 
Pinzke et al. [30], from which we derive the surface brightness 
of the clusters we analyze). r2oo = 1-9 Mpc (1.4°) for Perseus 
and 0.96 Mpc (2.8°) for Fornax [29]. The energy spectrum of 
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the model, in the energies of interests here (above 10 GeV), is 
a power-law with a slope of -2.25. 

Since the emission region is extended in the sky, we first 
divide the FOV into a grid of pixels each 0.2 degrees wide, 
and then we define the region of interest (ROI), constituted 
by all the pixels within an angle 9 max from the centre of the 
camera. We consider 15 values of energy threshold in log- 
arithmic steps from 50 GeV to 50 TeV. With the theoretical 
gamma-ray emission and the instrument response, we are able 
to compute the predicted number of background (A^° FF ) and sig- 
nal events (Nj) above each in each bin of the ROI sepa- 
rately, and then we integrate over the entire ROI. The model 
of Pinzke et al. [30] predicts a rather large gamma-ray flux 
for Perseus (O cr (> 100 GeV) = 2.04 x lO^'crrrV 1 ), the 
largest among the galaxy clusters, and a smaller one for Fornax 
(<D CR (> 100 GeV) = 1.5 x lO^cnrV 1 ). Above the different 
energy thresholds we determine how many hours CTA will 
need to detect the sources. We perform the calculation for the 
two CTA array B and E and for different ROI. We repeat the 
procedure 10 times for each energy threshold and average the 
results, in order to quantify the statistical fluctuations occurring 
when the number of events (both Nj and A^° FF ) are generated. 
The results are shown in Fig. 1.5. 
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Figure 1.5: Integration time required to have a 5<x detection (see text for 
details) of gamma-rays from CR-induced gamma-rays only according to the 
model of Pinzke et al. [30] for Perseus (lower curves, <I>cr(> 100 GeV) = 
2.04 X 10 _I1 crrr 2 s -1 ) and Fornax (upper curves, <£ C r(> 100 GeV) = 1.5 X 
10~ 13 cirr 2 s~'). The integration time is shown as a function of the energy 
threshold, and over different ROI, for the case of array B (dashed lines) and 
array E (solid lines). The shaded regions indicate the l<x standard deviation 
among 10 different simulations. 

If one assumes the CR-induced gamma-ray model by Pinzke 
et al. [30], CTA will detect such radiation from Perseus already 
in about 100 h, a fact which will constitute an extraordinary 
scientific result by itself 3 . The discovery could indeed be po- 



3 We underline that the upper limits obtained by the MAGIC experiment 



tentially close, opening up a completely new observation win- 
dow on the Universe. We underline that there is an absolute 
lower limit for gamma-rays in the hadronic scenario for clus- 
ters with an observed radio halo: a stationary distribution of CR 
electrons loses all its energy to synchrotron radiation for strong 
magnetic fields, as those in the radio halo, and therefore the ra- 
tio of gamma-ray to synchrotron flux becomes independent of 
the spatial distribution of the CRs and the thermal gas. For the 
Perseus cluster this lower limit is roughly a factor 3-4 from 
the gamma-ray flux predicted by the CR model [see Fig 3 in 
Ref. 47], hence CTA would, in the worst case scenario, require 
about 1 , 000 hours of observation to completely rule out the 
hadronic models. Such large observation times can in princi- 
ple be achieved either by, e.g., multi-annual observational cam- 
paigns. On the other hand, a non detection with CTA in a few 
hundred hours would seriously constrain the model and thus 
pose interesting challenges on the galaxy cluster physics. The 
situation is more pessimistic for our model of Fornax, which is 
out of reach for CTA. 

We see that the exact value of the integration time depends 
on the energy threshold chosen for the analysis. The reason for 
this is the tradeoff between the gamma-ray efficiency at differ- 
ent energies (the effective area), the source intrinsic spectrum 
and the chosen ROI. Roughly 90% of the CR-induced emission 
is expected within about 0. 1 r^o for Perseus, which corresponds 
to roughly 0.2°. We checked that integrating larger ROI, more 
background than signal is included in the analysis, thus dete- 
riorating the significance of the detection. This suggests that 
in realistic cases, the best ROI should be optimized. Finally, 
we also see that the prospects of detection are similar for both 
considered arrays, B and E. 

Gamma-ray emission from Dark Matter annihilation 

The gamma-ray brightness due to DM annihilations from a 
particular viewing angle in the sky is proportional to the DM 
density squared integrated along the line of sight, as shown in 
Eq. (1.1). In the case of galaxy clusters, the contribution of the 
smooth DM halo is boosted by the presence of DM subhalos. 
Recent A^-body simulations of Milky Way-like halos [5, 6, 7] 
found that the contribution of subhalos is small in the centre of 
the halo, due to dynamical friction and tidal effects that disrupt 
the subhalos. However, already at distances of 0.01 - 0.05 r2oo, 
subhalos become the dominant component. The real value of 
the boost factor from subhalos is unknown and the theoretical 
estimates depend on different assumptions and different meth- 
ods used in the calculations. Pinzke et al. [30] estimated a 
Bp = 580 and 910 for Fornax and Perseus respectively (for a 
minimal halo mass of 10~ 6 M Q ) , while other authors gave Bp 
from few tens [29] up to several thousands [31]. 

We refer again to the results of Pinzke et al. [30] where the 
authors assumed a double power-law to describe the luminosity 



on Perseus [47, Fig. 3] already constrain by about 20% the model predictions 
(the same used here), impling that the maximum CR acceleration efficiency is 
lower than 50% or, alternatively, the presence of non-neligible CR transport 
phenomena. 
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of subhalos as a function of the projected distance from the cen- 
tre of the halo, a behavior derived by analyzing the sub-halos in 
the Aquarius A^-body simulation. They also found the projected 
surface brightness to be largely independent of the initial pro- 
file of the smooth DM halo. As a result, the DM profile is very 
fiat since the emission decreases approximately only 10% at 
a distance of 1.5 - 2.0 degrees from the centre, depending on 
the cluster [Fig. 1.7 and Ref. 29]. For the case of Perseus and 
Fornax, we used the results of Fig. 10 of Ref. [30], assuming 
a telescope angular resolution of 0. 1 degree, which is a good 
approximation for CTA, despite the fact that the exact value 
depends on the array, the energy and the position in the FOV. 
We underline that in the case of galaxy cluster, the contribution 
from substructure strongly shapes the region of emission, basi- 
cally moving from a point-like source (in case no substructure 
are considered), to an extended source. Given the fact that the 
analysis differ in the two cases, the contribution from substruc- 
ture cannot be considered as a simple multiplicative factor in 
the intrinsic expected flux with respect to point-like case. 

Hereafter we consider the Fornax cluster, which has the 
largest expected DM-induced photon yield. The intrinsic flux 
is taken from [30, Table 2] and includes an intrinsic boost fac- 
tor from subhalos of Bp = 580, summing up to a total flux 
of <D DM (> 100 GeV) = 3.6 x lO^cnrV 1 . Additional intrin- 
sic boost factor may come from either other contributions from 
subhalos not accounted in this model, by mechanisms like the 
Sommerfeld enhancement discussed above, or by the effect of 
contraction processes due to baryonic condensation [110, 111]. 

To compute the CTA prospects of detection, we consider only 
the case of DM annihilating into bb (spectral shape obtained 
from Ref. [88]), while other channels like t + t~ or yfyT may be 
more constraining, depending on the energy (see Fig. 1.2). We 
take the reference thermal cross-section 3 x 10 _26 cm 3 s _1 and we 
scan DM particle mass m x between 50 GeV and 4 TeV. We op- 
timized the upper limit calculation as described in Ref. [38], by 
optimizing the energy threshold above which the upper limit is 
estimated. In addition, we consider the possibility of extending 
the size of the ROI up to a max of 2 degrees, to encompass the 
full radial extension of the source. Fig. 1.6 show the results. In 
100 h observation, the lack of detection would place exclusion 
limits at the level of 10~ 2 cm 3 s _1 . 

We also studied the effect of integrating over larger and larger 
regions: despite the increased numbers of background events, 
the signal yield is also larger and, in the case of Fornax, we gain 
more in integrating up to 9 max — 1° than 0.5°, while integrating 
over larger regions leads to a worse sensitivity. 

Distinguishing the dark matter signal from other gamma-ray 
contributions 

In the previous sections we have considered separately the 
contributions of CR and DM to the total gamma-ray photon 
yield. This is an unrealistic situation: galaxy clusters are, in 
fact, complex objects where gamma-rays may be due to dif- 
ferent contributions possibly of different spatial origin: by col- 
lisions of accelerated CRs, by DM annihilations and by fore- 
ground or embedded astrophysical sources. 
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Figure 1.6: Prospect of detection of DM-induced signal from Fornax for a DM 
particle annihilating into bb and 100 h integration time. The reference model is 
taken from Ref. [30] with subhalo boost factor Bp = 580. The shaded regions 
indicate the l<x standard deviation among 10 different simulations. 
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Figure 1.7: The surface brightness (above 1 GeV) of the gamma-ray emission 
from the Fornax cluste from CRs (red), DM (blue) and the sum of the two 
contributions (black). The DM emission is calculated from the K' benchmark 
model of [112] which has mass of 570 GeV and a velocity-averaged cross- 
section of 4.4 X 10~ 26 cm 3 s _1 . Adapted from Ref. [30]. 

Fortunately, gamma-rays of different origin typically have 
different spectral shapes, with the DM-induced emission char- 
acterized by the peculiar cut-off at E = m x and other remark- 
able spectral features [113, 1 12], in contrast to the plain spectral 
shapes (typically power-laws within the energy range of interest 
here) of the emission due to CRs, of the central galaxy or any 
astrophysical objects in the cluster. In the case a VHE emis- 
sion is detected from a cluster, this fact may be used as a probe 
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to discriminate between the components. However, we remark 
that in order to significantly discriminate the two sources one 
would need a quite significant detection over the CR-signal, 
which is often not supported by theoretical predictions for most 
galaxy clusters. 

A distinct approach could be based on the different spa- 
tial extensions of the various contributions of VHE gamma-ray 
photons from galaxy clusters. The possible individual galax- 
ies emitting within the cluster are typically seen as point-like 
sources, and thus one may exclude them from the FOV for CR 
and DM searches. Moreover, from the fact that CR-induced 
radiation is more concentrated than that induced by DM, one 
can optimize the ROI to select only those where the emission is 
DM dominated. In Fig. 1.7, we show the expected brightness 
profile for CR and DM photons for the Fornax cluster. One 
can see that up to 6 = 0.4° the emission is dominated by CR- 
induced photons, whereas this exact value is cluster-dependent 
and model-dependent, and in particular the possible intrinsic 
boost-factor in the DM signal can affect this. In this exam- 
ple, above 9 = 0.4°, the CR-signal fades more rapidly than 
the DM one. Then, in principle, by considering a ROI with 
a Omin = 0.4°, one could be able to isolate the DM signal. The 
maximum integration angle 8 max should be optimized accord- 
ing to the specific cluster and emission profile to maximize the 
sensitivity, as discussed above. Unfortunately, at the moment 
of writing this report, we did not have sufficient coverage in the 
MC of extended sources to perform such a study, and we are 
limited to a qualitative discussion. We mention that the "ge- 
ometrical" discrimination makes sense only if the DM signal 
is sufficiently large, otherwise different observational strategies 
could be more constraining. Finally, we stress again that with a 
large FOV (at least above 1 TeV) that has a near constant sen- 
sitivity over several degrees will allow CTA to study extended 
high energy gamma-ray sources in detail for the first time, with 
possibly revolutionary consequences for the IACT technique. 

1.3. Observations of the Galactic Halo and Centre 

The Galactic Centre (GC) is a long-discussed target for indi- 
rect DM searches with Cherenkov telescopes [114]. The den- 
sity of the DM halo should be highest in the very centre of the 
Milky Way, giving rise to a gamma-ray flux from annihilation 
of DM particles. On the one hand, this view is strengthened by 
the results of recent A^-body simulations of CDM halos [115] 
suggesting that, for an observer within the Milky Way, the an- 
nihilation signal from DM is not primarily due to small sub- 
halos, but is dominated by the radiation produced by diffuse 
DM in the main halo. On the other hand, searches close to 
the GC are made difficult by the presence of the Galactic Cen- 
tre source HESS J1745-290 [116, 117] and of diffuse emis- 
sion from the Galactic plane [118]. Both emissions can be 
plausibly explained by astrophysical emission processes: HESS 
J1745-290 is thought to be related to the Black Hole Sgr A* 
or the pulsar wind nebula G 359.95-0.04 [119], and the diffuse 
emission is well described as arising from hadronic cosmic rays 
interacting in giant molecular clouds. In both cases, the mea- 
sured energy spectra do not fit DM model spectra [120] and thus 



make a dominant contribution from DM annihilation or decay 
unlikely. 

In this situation, DM searches should better target regions 
which are outside the Galactic plane and hence not polluted 
by astrophysical gamma-ray emission, but which are still close 
enough to the GC to exhibit a sizable gamma-ray flux from DM 
annihilation in the Milky Way halo [121]. Given the angular 
resolution of Cherenkov telescopes and the scale height of the 
diffuse emission from the Galactic plane these criteria are ful- 
filled for an angular distance of about 0.3° from the GC. This 
angular scale translates into a distance of 45 pc from GC when 
using 8.5 kpc as the galactocentric distance. The radial DM 
density profiles obtained in N-body simulations of Milky Way 
sized galaxies, like Aquarius [115] and Via Lactea II [122], 
can be described by Einasto and NFW parameterizations, re- 
spectively. These parameterizations differ substantially when 
extrapolating to the very centre of the Milky Way halo since 
the NFW profile is much more strongly peaked. At distances 
greater than about lOpc, the difference is, however, just a fac- 
tor of 2 which implies that a search at angular scales of > 0.3° 
will not be hampered by the imprecise knowledge of the DM 
density profile at small scales. 

A search for a DM annihilation signal from the halo at angu- 
lar distances between 0.3° and 1.0° from the GC has recently 
been performed using 1 12h of H.E.S.S. data [42]. For WIMP 
masses well above the H.E.S.S. energy threshold of 100 GeV 
this analysis provides the currently most constraining limits 
on the velocity averaged annihilation cross section (cr ann v) of 
WIMPs (for IACTs) at the level of few 10~ 25 cmV 1 . Towards 
lower WIMP masses, observations of dwarf galaxies with the 
Fermi-LAT satellite yield even better limits [49] demonstrating 
how both observations of dwarf galaxies and of the extended 
GC region allow to jointly constrain the parameter space. 

Simulations and Assumptions 

The prospects of a search for DM annihilation photons from 
the Milky Way halo with CTA depend on (/) the performance 
of the southern CTA array, (ii) the applied analysis and back- 
ground rejection techniques, and (Hi) the details of the DM dis- 
tribution and WIMP annihilation. At low energies, the sensitiv- 
ity of IACTs is limited by the presence of hadron and electron 
showers which arrive isotropically and which can only be dis- 
tinguished from photons on a statistical basis. The basic strat- 
egy for the halo analysis is therefore to compare the fluxes of 
gamma-like events from a signal region (with solid angle AO s ) 
and a background region (solid angle AQ/,) and to search for 
DM features in the background-subtracted energy spectra. The 
signal region can be chosen such that it has the same instrumen- 
tal acceptance as the background region, but is located closer to 
the GC and features therefore a higher DM annihilation flux. 
For the purpose of this section, we rewrite Eq. (1.1) in terms of 
differential DM photon rate expected from the signal or back- 
ground regions (s, b respectively), given by: 

dR {o- wa v)dN 7 C 

-7ds,b= T^TTi 2 m)A(n,E)dQ, (1.3) 
dE %nm 2 x dE y J Af!sb 
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where dN y /dE y is the photon spectrum generated in the anni- 
hilation of a WIMP of mass m x , and A(fi, E) are the CTA ef- 
fective areas for photons, which depend on the position of the 
region within the FOV (£2), the energy E and further parameters 
(like the zenith angle of the observations). is the line-of- 
sight integral over the squared DM density p(r) (cf. Eq. 1.2). 
Since the DM density depends only on the distance to the GC 
r the line-of-sight integral and the astrophysical factor are only 
a function of the angular distance iff from the GC. Assuming 
that the signal and background region differ only with respect to 
their DM annihilation flux and their relative size a - AQ s /AOb, 
the rate of excess photon events R s - aR\, is given by 



Yl n E ^i \ [ 
l Jo dE y 



J(jff)A(Q,, E)dQ - 
J An, 



J AC. 



J(ifr)A(£l, E)dQ 
(1.4) 

Clearly, the rate vanishes when the astrophysical factors of 
the signal and the background regions are identical which im- 
plies that in the case of an isothermal DM density profile, a halo 
analysis with signal and background region chosen too close to 
the GC will not allow the placement of limits on (<x ann v). 




Figure 1.8: Illustration of the Ring Method for constructing signal and back- 
ground regions within one FOV of the CTA candidate arrays. The red star 
denotes the position of the GC in galactic coordinates; the blue star marks the 
pointing position of the CTA array which is shifted by an amount b in latitude 
from the GC. The annulus with inner and outer radii r\ and ri around the ob- 
servation position defines regions of equal acceptance. The signal region (blue, 
close to the GC) is constructed as intersection of the annulus and a circle around 
the GC with radius A cut . The remaining regions on the annulus (red) are used 
as background region. Regions within ±0.3° of the galactic plane (yellow) are 
neither part of the signal nor of the background region. 
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Table 1.2: Optimized values of the parameters used in the application of the 
Ring Method for the candidate arrays E and B. See Fig. 1.8 for a description of 
the parameters. 



Given an observation time T, Eq. 1.4 can be used to esti- 
mate the number of excess photons for a particular realization 
of CTA and a DM model defining (cr ann v), dN y fdE y and J(iff). 
Equivalently, one can place a limit on (cr ann y) given an upper 
limit on the number of excess photon events. Simulations of 
the candidate arrays E and B at a zenith angle of 20° were used 
to infer the effective area for diffuse photons and the residual 
rate of protons anywhere in the FOV. Both arrays feature large- 
size telescopes and are therefore suitable for studies in the low- 
energy domain. The available observation time was set to 100 h, 
which is about 10 % of the total observation time per year. 

Two different ways of defining signal and background re- 
gions were employed and compared, namely the so-called Ring 
Method and the On-Off Method. For the Ring Method, the can- 
didate arrays E or B were assumed to observe the GC region at 
Galactic longitude / = and Galactic altitude b, and signal and 
background regions were placed in the same FOV as illustrated 
in Fig. 1.8. An annulus with inner radius r\ and outer radius 
ri around the observation position was constructed and divided 
into signal and background region such that the signal region is 
closer to the GC and has therefore a larger astrophysical factor. 
The separation of signal and background region is achieved by 
a circle with radius A cut around the GC whose intersection with 
the annulus defines the signal region. All other regions on the 
ring were considered as background region. The values of the 
four parameters b, r\ , ri and A cut were optimized such that the 
attained significance of a DM signal per square root time was 
maximized. The maximization was carried out for a wide range 
of WIMP masses but the dependence on the actual WIMP mass 
was found to be fairly weak. The derived values for both can- 
didate arrays are listed in Tab. 1.2. Judging from present IACT 
observations, we do not expect strong diffuse gamma-ray emis- 
sion to extend outside the +0.3° box used to mask the galactic 
disc. New point-like or slightly extended sources will be ex- 
cluded, making the On and Off region smaller. In addition, the 
approach is only sensitive to gradients in the diffuse gamma-ray 
emission, whereas the charged particle background is isotropic. 
In the optimization process an Einasto profile was assumed for 
the DM signal, but the optimal values are only weakly depen- 
dent on the assumed profile in the region beyond 0.3 degrees 
from the Galactic plane. 

The usage of the annulus implies the same acceptance for 
signal and background region since the acceptance is, to good 
approximation, only a function of the distance to the observa- 
tion position. Placing both signal and background regions in the 
same FOV implies that both regions will be affected by time- 
dependent effects in a similar way. A disadvantage is, however, 
that the angular distance between the signal and background re- 
gion is only of order of the FOV diameter, reducing the contrast 
in Eq. 1 .4 considerably. This contrast was increased in the On- 
Off Method where data-taking with an offset of typically 30' in 
Right Ascension was assumed. In this mode, the telescopes first 
track for half an hour the same observation position as in the 
Ring Method which defines the signal region. The telescopes 
then slew back and follow the same path on the sky for another 
30min. The second pointing has the same acceptance as the 
first one since the same azimuth and zenith angles are covered 
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but generates a background region with much increased angu- 
lar distance to the GC. In the On-Off Method, the observation 
time was 50 h for the signal and 50 h for the background region 
giving again a total observation time of 100 h. Regardless of 
whether the Ring Method or the On-Off Method was used, all 
areas with \b\ < 0.3° were excluded from signal and background 
regions to avoid pollution from astrophysical gamma-rays. 
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Table 1.3: Astrophysical factor for the signal region and size of the integration 
region (Af2) for Ring and On-Off method. In the case of the On-Off Method, 
An was chosen as the entire FOV of the candidate array which introduces a 
dependence on the assumed WIMP mass since the effective FOV grows with 
photon energy. The table gives values for a WIMP mass of 0.1, 1, and 10 TeV. 

The astrophysical factor (Eq. 1 .2) was taken from the Aquar- 
ius Simulation [115] which had been corrected for the pres- 
ence of subhalos below the resolution limit of the simula- 
tion. The line-of-sight integral assumes a value of 40.3 x 
10 24 GeV 2 cnr 5 sr 1 at ¥ = 1°. Table 1.3 lists the astrophysical 
factors of the signal regions which were defined in the Ring and 
On-Off Method, respectively. In case of the On-Off Method, the 
signal region was defined as the total effective FOV of the On- 
pointing which introduces a dependence on the WIMP mass 
since the FOV grows with photon energy. For the WIMP anni- 
hilation spectrum dN 7 /dE y several different choices were con- 
sidered. The generic Tasitsiomi spectrum [60] is appropriate for 
a dominant annihilation into quark-antiquark pairs with subse- 
quent hadronization into n° particles and was used in the opti- 
mization of the parameters of the Ring Method. Other spectra 
were explored by considering bb, t + t~ and yU + //~ final states 
[88]. 

Discussion 

The two plots in Fig. 1.9 show the upper limits for WIMP 
masses between 0.1 TeV and 10 TeV, translated from the sensi- 
tivity using here the method of Feldman and Cousins [59]. Each 
curve corresponds to one set of assumptions. It is evident that 
the most constraining limits can be derived for masses of about 
0.5 TeV which is a factor of 2 improvement compared to cur- 
rent IACT arrays like H.E.S.S. reaching best sensitivity around 
1 TeV. This is a direct consequence of the lower threshold and 
superior stereoscopic background rejection of the CTA candi- 
date arrays. Typical limits are around few 10 _26 cm 3 s _I which 
is a factor of 10 improvement compared to current IACTs. The 
comparison of array E (blue) and B (same line style but red) 
shows that the limits for array B are always better, which can 
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Figure 1.9: CTA sensitivities on the velocity averaged annihilation cross- 
section as a function of the WIMP mass. Shown are curves for the candidate 
arrays E (blue) and B (red). Top: Comparison of the Ring Method (solid lines) 
and On-Off Method for background subtraction. Annihilation as in Tasitsiomi 
and Olinto [60] was assumed. Bottom: Comparison of different WIMP spectra 
for the Ring Method. The solid line denotes the case of annihilation into bb; 

and t + t~ spectra are shown by the dotted and dashed lines, respectively. 
On both panels, the classical annihilation cross section for thermally produced 
WIMPs at 3 X 10 -26 cm 3 s~' is indicated by the black horizontal line. 



be understood from the fact that B comprises 5 large-size tele- 
scopes and array E only 4. The magnitude of this effect is, 
however, comparatively small (~ 20 %). Overall, CTA should 
be able to probe the parameter space below the velocity aver- 
aged annihilation cross-section for thermally produced DM of 
3 x 10~ 26 cm 3 s~ 1 for WIMP masses between several ten GeV 
and several TeV. 

The upper panel of Fig. 1.9 illustrates the impact of data- 
taking with the Ring Method and the On-Off Method for the 
case of a dominant annihilation into quark-antiquark pairs with 
subsequent 7r° creation [60] . The On-Off Method (dashed lines) 
is more sensitive than the Ring Method (dashed lines). One 
must keep in mind, however, that the On-Off Method spends 
50 % of the observation time far away from the GC which im- 
plies that this data set will be of limited use for studies of astro- 
physical sources. Another drawback of the On-Off Method is 
its susceptibility to systematic effects arising from variations of 
the data-taking conditions (electronics, atmosphere). In view of 
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this, the increased sensitivity for the DM halo analysis in parts 
of the parameter space will not probably suffice to motivate the 
acquisition of a larger data set in this mode. 

Compared with the choice of the CTA candidate array (B or 
E) and the analysis method (Ring Method or On-Off Methods), 
the WIMP annihilation spectrum has the strongest impact on 
the CTA sensitivity. The lower panel of Fig. 1 .9 shows for both 
candidate arrays and the Ring Method the limits obtained in the 
case of a dominant annihilation into bb pairs (solid), /i 4 pT (dot- 
ted) and t + t~ (dashed). The small photon yield from p + p~ final 
states implies limits that are a factor of about 5 worse than lim- 
its for dominant annihilation into t + t~. It is clear that the full 
potential of the halo analysis will be exploited by confronting 
individual DM models with their predicted WIMP annihilation 
spectra dN y /dE y with data. 

1.4. Anisotropies in the diffuse gamma-ray background 

Besides gamma-rays from individual resolved sources and 
Galactic foreground, another component of diffuse gamma-ray 
background radiation has been detected and proven to be nearly 
isotropic. This radiation dominantly originates from conven- 
tional unresolved point sources below the detection threshold, 
while another fraction might be generated by self-annihilating 
(or decaying) DM particles, which then could produce spe- 
cific signatures in the anisotropy power spectrum of the diffuse 
gamma-ray background [123, 124, 125, 126, 127]. The differ- 
ent hypotheses about the origin of the gamma-ray background 
may be distinguishable by accurately measuring its anisotropy 
power spectrum. 

Compared to the current generation of IACTs, CTA will have 
improved capabilities to measure anisotropies in the diffuse 
gamma-ray background, based upon a better angular resolution 
(determined by the point-spread function, PSF), an increased 
size of the FOV, and a higher background rejection efficiency. 
In the following, we discuss the effects of different assumptions 
on the background level and the anisotropy spectrum on the re- 
construction of the power spectrum for the current generation 
of IACTs, and address the improvement obtainable with CTA. 
Finally, we make predictions for the discrimination between as- 
trophysical and dark matter induced anisotropy power spectra 
for CTA. 

Simulation 

In order to investigate the measured power spectrum and 
the impact of instrumental characteristics, a sample of event 
lists containing anisotropies generated with Monte-Carlo simu- 
lations was analyzed. The event lists were simulated by gener- 
ating skymaps with a given anisotropy power spectrum. In total, 
12 skymaps covering the size of the FOV and being in different 
celestial positions were created, with a power spectrum for a 
given multipole moment I defined as Q = l/(2{ + 1) Y} a (m\ 2 , 
m = -I, ...,{, where ae m denotes the coefficients of a (real- 
valued) spherical function decomposed into spherical harmon- 
ics. With {a tm ) = 0, C( reflects the width of the a lm distri- 
bution, which was assumed to be Gaussian. The simulations 
were made for different power spectra l(t + 1)Q ~ £ s , with 



s = 0.5, 1.0, 1.5, 2.0, 2.5. The pixel size of these skymaps was 
0.002°, corresponding to I = 9 x 10 4 (where ® t = 180°/£). The 
skymaps /(#, <f) were normalized in a way that the pixel with 
the smallest signal was assigned the value and the pixel with 
the largest signal was assigned 1. Anisotropy power spectra 
were then derived from the fluctuation maps I(§,(p)/(I), such 
that for a full signal the maximum allowed difference in each 
map equals 1 . Note that this difference can be smaller when an 
additional isotropic noise component is present. 

An event was simulated in three subsequent steps: First, the 
celestial position was randomly chosen within the FOV, and the 
event was classified to represent a signal- or isotropic noise- 
event, respectively. The decision for a signal event was based 
upon a normalized random number z: If z was smaller than 
the skymap value at the corresponding position, the event was 
considered a signal event. Otherwise, another event position 
was selected while reapplying the procedure. Subsequently, the 
event map was convolved with a PSF of 0.1°, which is similar to 
the resolution of current IACTs. The effect of a better angular 
resolution is discussed below. The event maps were simulated 
to contain 10 7 entries. Note here that this number, as selected 
for the toy model, does in general not reflect the actual number 
of expected physical signal events. Therefore, the following 
discussion is focussed more on a qualitative discussion of the 
criticalities of the calculation rather than on making quantitative 
predictions. 

To analyze an event list containing N &v events, a HEALPix 
skymap with N p i x pixels was accordingly filled, and analyzed 
using the HEALPix software package 4 . Therefore, the analyzed 
function is 

N ■ iVpix 

w (h)f(h) = -^Yx i -b i (h), (1-5) 

1 Vevent ._ j 

where x, denotes the number of events in pixel i, and fe,-(n) 
equals 1 inside pixel ; and outside. The function w(n) de- 
scribes the windowing function — in this case the FOV with 
Gaussian acceptance — and f(h) denotes the original signal 
function over the full sky. The windowing function was nor- 
malized such that the integral over the full sky equals An: 

j dQ.w(h)f(h) = 4n. (1.6) 

Note that this differs from other analyses of this type, where w 
is defined such that the maximum value is 1 . This difference in 
the normalization was done in order to keep a simple simula- 
tion code, and results should be equivalent. Final results were 
averaged over the corresponding 12 skymaps. 

The effect of the anisotropy spectrum and the residual back- 
ground level on the spectral reconstruction 

In Fig. 1.10, we show the mean value and the RMS of the 
Ce power spectra. The value Q represents the strength of 
anisotropies of the angular scale = 180°/^. Anisotropies 
smaller than the angular resolution (defined by the PSF) are 



4 http://healpix.jpl.nasa.gov/ 
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Figure 1.10: Measured power spectra /(/ + l)Cg/(2n) for different slopes s of 
the simulated input spectrum, compared to an isotropic background spectrum. 
Color-filled areas depict the RMS of the spectra. The size of the PSF is cxpsf = 
0.1°. For reference, the simulated spectra are shown as dashed lines. 



Figure 1.11: Influence of the signal fraction on the measured power spectrum. 
Shown are the reproduced spectra for a slope * = 0.5 for several ratios between 
signal and total events; the background events are distributed isotropically. In 
order to estimate the effect of the noise ratio, the best fit levels are shown as 
dotted lines. The width of the PSF is chosen as in Fig. 1.12, crpsF = 0.05°. 



smeared out. This effect is clearly visible for large ( > 1000, 
where the power spectra converge into the Poissonian noise of 
the isotropic background spectrum. The angular resolution as- 
sumed for the simulation shown in this figure has a width of 
o"psf = 0.1°. Furthermore, anisotropies with a size larger than 
the FOV are truncated at t ~ 100 due to the effect of the win- 
dowing function. The simulated FOV in Fig. 1.10 has a width 
of 2.5°, which is comparable with the FOV of current IACT ex- 
periments. For the toy model, Fig. 1.10 demonstrates that, for 
t ~ 100-1000, power spectra of different slopes are separable 
within the statistical errors and distinguishable from isotropic 
noise. CTA will have a smaller PSF as well as a larger FOV. 
This will make the signal vanish at larger I than in the exam- 
ple, and the windowing function will influence the spectrum to 
smaller I than in the figure. Therefore, we conclude that the 
FOV as well as the PSF, while important, will not be crucial for 
the investigation of anisotropies with CTA in the desired multi- 
pole range. 

In general, the measured flux will be composed of both signal 
and background events. The background is produced mainly by 
two separate processes: 

1. Events caused by cosmic rays (protons and electrons) 
which are misinterpreted as photon events. 

2. An isotropic component of the photon background radia- 
tion, which does not count as signal according to our defi- 
nition. 

The influence of isotropic background is demonstrated in 
Fig. 1.11, where the power spectrum for s = 0.5 is shown for 
different background levels. Here, the signal fraction is defined 
by /sig = N s ig/N eY , where A^ s i g denotes the number of signal 
events. The overall power is clearly reduced in case of fully 
isotropic background. From the figure, we see that when the 
signal fraction improves by a factor 5, the power spectrum is 
boosted by about two orders of magnitude. For this reason, we 
expect the ten-fold improved CTA sensitivity to mark the major 
difference with respect to the current generation of IACTs for 



such studies. 
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Figure 1.12: Comparison between measured anisotropy power spectra with (a) 
a pure astrophysical origin represented by CV = 10~ 5 (blue bands) and (b) 
with an additional DM component, i.e., 20% of the total flux, represented by 
C{ = 10~ 3 (red bands). The assumed observation time is 300 h. The three cases 
in each plot represent the hadronic background rates of 10 Hz, 1 Hz, and 0. 1 Hz. 

The theoretical expectations for the power spectra of the dif- 
fuse gamma-ray flux of both the astrophysical as well as the 
DM components are highly model dependent. Since the astro- 
physical component is dominated by the gamma-ray flux from 
unresolved point sources, expected with a constant Q (s = 2.0 
in our notation), we conservatively assume the slope of the 
DM component(s) to be similar. In this scenario, the differ- 
ence between the power spectra manifests in the normaliza- 
tion. For unresolved point sources, C^biazars = 10~ 5 , while for 
DM-induced anisotropies, considering the thermal annihilation 
cross-section 3 x 10~ 26 cm 3 s _1 , Q_dm = 10~ 3 is expected [see, 
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e.g., 127]. In our simulation, this was realized by distributing 
N = 4n/C( point sources over the full sky. While representing 
a non-physical model, this is a convenient way of producing a 
Poissonian anisotropy power spectrum which is a reasonable 
assumption for generic astrophysical and DM emitters. The 
normalization of the signal was set by extrapolating the spec- 
trum of the extragalactic gamma-ray background (EGB) [128] 
to E > 100 GeV. Note that the strength of the DM annihi- 
lation signal is strongly affected by the formation histories of 
DM halos and the distribution of DM subhalos. For example, 
Fig. 3 in [129] shows that the gamma-ray spectrum of DM 
annihilation could reach the measured gamma-ray background 
spectrum and therefore deliver a significant fraction of the mea- 
sured flux. Here, we investigate the cases that (a) the total EGB 
originates from astrophysical sources and (b) 20% of the EGB 
(optimistically) originates from DM annihilation. The isotropic 
hadronic component depends on analysis cuts and the quality of 
the gamma-hadron separation. In the following, three different 
background rates are assumed: 10 Hz, 1 Hz, and an optimistic 
0.1 Hz rate. We assume a CTA-like FOV with a radius of 5° and 
a CTA-like PSF with cr PSF = 0.05°. The results are shown in 
Fig. 1.12, where each band represents a sample of 20 realiza- 
tions. One can see in the figure that depending on the achieved 
background rate, in principle the two above mentioned scenar- 
ios (a) and (b) will be well distinguishable for CTA. 

2. Search of axion-like particles with CTA 

Axions were proposed in the 1970's as a by-product of the 
Peccei-Quinn solution of the strong-CP problem in QCD [130]. 
In addition, they are valid candidates to constitute a portion of, 
or perhaps the totality of, the non-baryonic CDM content pre- 
dicted to exist in the Universe. Another extremely interesting 
property of axions, or more generically, Axion-Like Particles 
(ALPs, for which — unlike axions — the mass m a and the cou- 
pling constant are not related to each other), is that they are 
expected to convert into photons (and vice versa) in the pres- 
ence of magnetic fields [131, 132]. The photon/ALP mixing is 
indeed the main signature used at present in ALP searches, such 
as those carried out by CAST [133, 134, 135, 136] or ADMX 
[137], but it could also have important implications for astro- 
nomical observations. For example, photon/ALP mixing could 
distort the spectra of gamma-ray sources, such as Active Galac- 
tic Nuclei (AGN) [138, 139, 140, 141] or galactic sources, in 
the TeV range [142]. 

The photon/ALP mixing effect for distant AGN was also 
evaluated by Sanchez-Conde et al. [143] under a consistent 
framework, where mixing takes place inside or near the gamma- 
ray emitter as well as in the intergalactic magnetic field (IGMF). 
A diagram that outlines this scenario is shown in Fig. 2.1. The 
artistic sketch shows the travel of a photon from the source to 
the Earth and the main physical cases that one could identify 5 . 



5 Note that this formalism neglects, however, the mixing that may hap- 
pen inside the Milky Way due to galactic magnetic fields. In the most ide- 
alistic/optimistic case, it would produce a photon flux enhancement at Earth 
of -3% [141]. 



From top to bottom: 1) the photon converts to an axion and 
back to photon in the IGMF, 2) the photon converts to an axion 
in the IGMF, 3) the photon converts to an axion at the source, 
which then does not interact with the EBL, therefore traveling 
unimpeded from the source to the Earth, 4) the photon travels 
unimpeded from the source to the Earth, 5) the photon converts 
to an axion at the source and back to photon in the IGMF, 6) the 
photon interacts an EBL photon resulting in pair production. It 
is clear that cases 2, 3, 6 corresponds to an attenuation of the 
intrinsic source flux, while cases 1,4,5 allow for a recovery of 
the intrinsic photon yield. 




Figure 2. 1 : Photon/ALP conversions (crooked lines) that can occur in the emis- 
sion from a cosmological source, y and a symbols represent gamma-ray pho- 
tons and ALPs respectively. This diagram collects the main physical scenarios 
that we might identify inside our formalism. Each of them are schematically 
represented by a line that goes from the source to the Earth. From [143]. 

The probability of a photon of energy E y to be converted into 
an ALP (and vice versa) can be written as [138]: 



Po = 



1 



1 + (E cri ,/E r ) 2 




(2.1) 



where s is the length of the domain where there is a roughly 
constant magnetic field B, and M the inverse of the coupling 
constant. Here we also defined a characteristic energy, E cr i t \ 
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or in more convenient units: 



E crit (GeV) = 



m 2 „ Mn 

[ieV 1 1 

0.4 B c 



(2.2) 



(2.3) 



where the subindices refer to dimensionless quantities: 
m^ e v = w/jueV, M u = M/10 11 GeV and B c = B/Gauss; 
m is the effective ALP mass m 2 = \m 2 - u> 2 pl \, with 

copi = 0.37 x lO^VeV -<Jn e /cm~ 3 the plasma frequency 
and n e the electron density. The most recent results from 
the CAST experiment [136] give a value of M u > 0.114 
for ALP mass m a < 0.02 eV. At present, the CAST bound 
is the most general and stringent limit in the range 10" 11 eV 
«: m a «: 10~ 2 eV. 

The main effect produced by photon/ALP mixing in the 
source is an attenuation in the total expected intensity of the 
source just above a critical energy E cr j t (see Fig. 2.1). As for 
the mixing in the IGMFs, despite the low magnetic field B, 
the photon/ALP conversion can take place due to the large dis- 
tances involved. In the model of Sanchez-Conde et al. [143], 
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it is assumed that the photon beam propagates over N domains 
of a given length. The modulus of the IGMF is the same in 
all of them, whereas its orientation changes randomly from one 
domain to the next, which in practice is also equivalent to a 
variation in the strength of the component of the magnetic field 
relevant to the photon/ALP mixing. 

In discussing photon/ALP conversion in IGMFs, it is also 
necessary to consider the important role of the Extragalac- 
tic Background Light (EBL), its main effect being an addi- 
tional attenuation of the photon flux (especially at energies 
above about 100 GeV). Recent gamma-ray observations al- 
ready pose substantial challenges to the conventional models 
that explain the observed source spectra in terms of EBL atten- 
uation [144, 145, 146, 147, 148]. 

Taken together, photon/ALP conversions in the IGMF can 
lead to an attenuation or an enhancement of the photon flux 
at Earth, depending on distance, magnetic fields and the EBL 
model considered. A flux enhancement is possible because 
ALPs travel unimpeded through the EBL, and a fraction of them 
can convert back into photons before reaching the observer. 
Note that the strength of the IGMFs is expected to be many or- 
ders of magnitude weaker (~nG) than that of the source and its 
surroundings (~G). Consequently, as described by Eq. (2.3), the 
energy at which photon/ALP conversion occurs in this case is 
many orders of magnitude larger than that at which conversion 
can occur in the source and its vicinity. Assuming a mid-value 
of B-0.1 nG, and M n =0.114 (CAST lower limit), the effect 
could be observationally detectable by IACTs only if the ALP 
mass is of the order of 10~ 10 eV, i.e. we need ultra-light ALPs. 

In order to quantitatively study the effect of photon-axion 
conversion over the cosmological distances of AGN, we con- 
sider the total photon intensity. It becomes then useful to define 
the axion boost factor as the difference between the predicted 
arriving photon intensity without including ALPs and that ob- 
tained when including the photon/ALP conversions. Qualita- 
tively speaking, it is found that the more attenuating the EBL 
model considered, the more relevant the effect of photon/ALP 
conversions in the IGMF (since any ALP to photon reconver- 
sion might substantially enhance the intensity arriving at Earth). 
Furthermore, higher B values do not necessarily translate into 
higher photon flux enhancements. There is always a B value 
that maximizes the axion boost factors; this value is sensitive 
to the source distance, the considered energy and the adopted 
EBL model (see Ref. [143] for a more detailed discussion). 

There could be indeed different approaches from the observa- 
tional point of view, although all of them will be probably based 
on the search and analysis of a systematic residual after apply- 
ing the best-fit (conventional) model to the AGN data. For ex- 
ample, Ref. [143] predicts the existence of a universal feature in 
the spectrum of the sources due to the intergalactic mixing, that 
is completely independent on the sources themselves and only 
depends on the ALP and IGMF properties. This feature should 
be present at the same critical energy E crit for all sources, and 
would show up in the spectra as a drop in the flux — whenever 
E crit is in the range where the EBL effect is negligible — or 
even as a sudden flux increase, if the EBL absorption is strong 
for E = Ecrit. 



Test case for CTA: PKS 1222+21 

We have taken as a test source the flat spectrum radio quasar 
4C +21.35 (PKS 1222+21), at redshift z = 0.432, which was 
detected by MAGIC above 70 GeV [149] in June 2010, dur- 
ing a target of opportunity observation triggered by the high 
state of the source in the Fermi-LAT energy band. This source 
is the second most distant object detected by ground-based 
gamma-ray telescopes, and hence an ideal candidate for the 
study of propagation effects. The observed energy spectrum of 
4C +21.35 during the 0.5 hour flare recorded by MAGIC was 
well described by a power law of index F = 3.75 ± 0.27 sm ± 
0.2 syst . The intrinsic spectrum, assuming the EBL model of 
Dommguez et al. [150] was estimated to be a power law of in- 
dex F = 2.72 ± 0.34, which extrapolated down to an energy 
of about 5 GeV, connects smoothly with the harder spectrum 
(r = 1.95 ± 0.21) measured by Fermi-LAT between 0.2 and 
2 GeV in a 2.5 h period encompassing the MAGIC observation. 
It must be noted that longer-term Fermi-LAT observations of 
the source in various states of activity show a break in the spec- 
trum between 1 and 3 GeV, with a spectral index after the break 
(and up to 50 GeV) ranging between 2.4 and 2.8 [151]. 

We have simulated CTA observations of 4C +21 .35 assuming 
an intrinsic unbroken power-law spectrum, in the relevant en- 
ergy range, like the one determined by MAGIC for 4C +21.35 
during the flare, i.e. dN/dE = K x[£/(0.2 TeV)] -272 . Keep- 
ing the spectral shape unchanged, we have tried different abso- 
lute flux normalizations, taking as a reference the flux observed 
by MAGIC, K = 1.78 x 10~ 5 nr 2 s" 1 TeV" 1 . We have also 
tested different observation times: the actual duration of the 
VHE flare observed by MAGIC is unknown, since the observa- 
tion was interrupted while the flare was still going on, but the 
flares observed by Fermi-LAT above 100 MeV show rise and 
decay time scales of the order of a day [151], so it is reasonable 
to expect that the source may stay several hours in flux states as 
high as that observed by MAGIC. For the detector simulation 
we have used the CTA candidate array E. The EBL model in 
Ref. [150] has been used to account for the effect of the EBL, 
and the conversion of photons into ALPs and vice versa has 
been simulated following the formalism detailed in Ref. [143] 
as outlined above. Only conversions in the IGMF have been 
considered (in this case, mixing in the source typically leads 
to only a few percent of flux attenuation, so we neglected it in 
order to avoid extra uncertainties). 

We have assumed the same parameters for the IGMF as those 
in the fiducial model in [143]: 0.1 nG is the (constant) modulus 
of the IGMF 6 , which is assumed to have fixed orientation within 
domains of size 1 Mpc. The orientation of the IGMF varies 
randomly from one domain to the next. The ALP parameters, 
mass and coupling constant, enter via E cr i t , below which the 
conversion probability is negligible. We have scanned E crit in 
the range 0. 1 to 1 TeV in steps of 0. 1 TeV 



6 Note that this is just one order of magnitude below the current upper limits 
[152], and lower limits lie many orders of magnitude below [153, 154]. With a 
IGMF of 0.01 nG or smaller the effect of the photon/ALP mixing would prob- 
ably be too weak to be observed with CTA. 
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Figure 2.2: Simulation of a 5 h CTA observation of a 4C +21.35 flare 5 times 
more intense than the one recorded by MAGIC [149]. In black, energy bins 
used for the fit (those with a signal exceeding three times the RMS of the back- 
ground, and a minimum of 10 excess events). Excluded points are displayed 
in grey. The estimated intrinsic differential energy spectrum (after correcting 
for the EBL effect) shows a boost at high energies due to photon/ALP mix- 
ing. The IGMF strength is assumed to be 0.1 nG, and ALP parameters result in 
En* = 200 GeV. 
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Figure 2.3: Same as in Fig. 2.2, but with E„-, f = 1 TeV. Note that in scenarios 
like this, where E cn - ( is within the energy range in which the EBL absorption is 
already large, the boost in the flux shows up as a sudden rise (smeared out by 
the spectral resolution of the instrument) which would even allow to determine 
E cr ,-, accurately. 



Using the performance parameters of array E, we obtain the 
expected gamma-ray and cosmic -ray background rates in bins 
of estimated energy, and from them the reconstructed differ- 
ential energy spectrum. After this, we correct the observed 
spectrum by the energy-dependent attenuation factors expected 
from the EBL in order to get an estimate of the intrinsic source 
spectrum. Each simulated spectrum is fitted to a power-law 
with variable index of the form dN/dE oc E'"^ '°s( £ /°- 1TeV ), 
in which we constrain the /3 parameter so that the spectrum can- 
not become harder with increasing energy (such behavior is not 
expected from emission models in this energy range). Only en- 
ergy bins with a signal exceeding three times the RMS of the 
background, and a minimum of 10 excess events, are consid- 



ered in the fit. 

In the absence of any significant photon/ALP mixing, the 
resulting fits will all match the spectral points within the ex- 
perimental uncertainties, resulting in good x 1 values. But, as 
shown in Ref. [143], certain combinations of ALP parameters 
and values of the IGMF may result in significant modifications 
of the observed VHE spectra. The most striking feature is a 
boost of the expected flux at high energies, which is particularly 
prominent in the estimated intrinsic (i.e. EBL-de-absorbed) 
spectrum. Such a feature may result in a low value of the x 2 - 
probability of the spectral fit. In Figs. 2.2 and 2.3 we show two 
such cases, in which the observed spectra, after de-absorption 
of the EBL effect, show a clear hardening of the spectral in- 
dex. The effect is particularly striking in the cases in which the 
EBL absorption at E = E cn -, is already strong (e.g. Fig. 2.3), be- 
cause then the boost sets in very fast, resulting in dN/dE rising 
with energy at around E tT „. The rise is actually very sharp, but 
it is smoothed by the energy resolution of the instrument. An 
improvement in the energy resolution would increase the sig- 
nificance of the feature and improve the determination of E crit . 
In contrast, if E c „, is in the range in which the EBL absorp- 
tion is small or negligible (Fig. 2.2), the feature at E cn -, would 
just be a flux drop of at most =s 30% [143], also washed out 
by the instrumental energy resolution. In those cases, though 
a high-energy boost may still be clearly detected, it would be 
hard to determine the exact value of E fr „. This is because, in 
the formalism described in Ref. [143], similar ALP boost fac- 
tors are always achieved at energies E > E cr i t , independently of 
the particular value of E„„ in each case. 
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Figure 2.4: Median of the % 2 -probabilities of the fits to the de-absorbed differ- 
ential energy spectra of 4C +21.35 measured by CTA, assuming photon/ALP 
mixing, for different values of E c „ ( . We simulated observations of flares of 
two different durations: 0.5 and 5 hours, and with intensities equal to 1 and 5 
times that of the flare reported in [149]. The dashed horizontal line marks the 
probability that corresponds to 5 standard deviations. 



Prospects 

For each of the E cr n values scanned, we have performed 10 3 
simulations of a CTA observation, all with the same source flux 
and observation time. We consider that a given value of E cr u 
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is within the reach of CTA whenever the median of the x 2 - 
probability distribution is below 2.9 x 1(T 7 , which corresponds 
to 5 standard deviations. In Fig. 2.4 we show the median of the 
X 2 probability versus E crit , for two different assumptions on the 
source flux and two different observation times. The range of 
E C rit which can be probed with CTA for the different scenarios 
is the one for which the curves in Fig. 2.4 are below the dashed 
horizontal line. As expected, the range becomes larger as we 
increase the observation time and/or the flux of the source. A 
0.5 h duration flare like the one reported in [149] would not 
be enough for CTA to detect a significant effect in any of the 
tested ALP scenarios, i.e. the solid black line never goes below 
the dashed line for any value of E cnf . A flare of similar inten- 
sity, but lasting 5 hours (green line) would already be enough to 
see the boost due to ALPs for those scenarios with E cn , < 500 
GeV. In Fig. 2.4 we can also see that for a hypothetical flare 
with an intensity 5 times larger, lasting 5 hours, the accessible 
range of E crit would extend up to 1.3 TeV. 

3. High Energy Violation of Lorentz Invariance 

Lorentz Invariance (LI) lies at the heart of all of modern 
physics, in particular the unification of space and time through 
the principle of Special Relativity. Space-time was elegantly 
promoted to be a dynamic entity in the covariant classical the- 
ory of gravity namely General Relativity (GR) which has been 
rigorously tested on astronomical scales and underlies the math- 
ematical description of cosmology. Similarly quantum me- 
chanics has been successfully married with Special Relativ- 
ity to yield the quantum theory of fields which underlies the 
very successful Standard Model of leptons and quarks and the 
gauged electromagnetic, weak and strong forces. However it 
has proved considerably more difficult to unify gravity with the 
other forces, since GR is fundamentally non-renormalisable. A 
fully quantum theory of gravity (QG) is still beyond our grasp 
although there has been significant progress towards this goal in 
various approaches such as superstring theory and loop quan- 
tum gravity [155]. QG should describe dynamics at the Planck 
energy £ P i = y/hc/G N - 1.22 x 10 19 GeV or equivalently the 
Planck length l PX = ^HG N /c 3 1.62 x 10 cm, where gravi- 
tational effects should become as strong as the other forces and 
the notion of space-time is likely to need revision. This has 
opened up the possibility that LI may be violated by QG effects 
although, lacking a fully dynamical theory, the expectation is 
generic rather than definite. For example quantum fluctuations 
may produce 'space-time foam' at the Planck scale resulting in 
a non-trivial refractive index and anomalous dispersion of light 
in vacuo i.e. an energy dependence of the speed of light. Hence 
over the past decade, there has been tremendous interest in test- 
ing LI at high energies as part of what has come to be called 
'quantum gravity phenomenology' [156, 157, 158]. 

Possible energy dependence of the speed of light in the vac- 
uum has been predicted, in the framework of several theories 
dealing with quantum gravity models and effective field theory 
models [159]. The seminal paper by Amelino-Camelia et al. 
[160] proposed that this can be parameterised by a Taylor ex- 



pansion of the usual dispersion relation: 

c 2 p 2 = E 2 [l ± ti(E/E n ) ± &(EIE n ? + ■■■], (3.1) 

where the value of the co-efficients % a would be specified by the 
theory of quantum gravity (and may well turn out to be zero) . 
For example there are specific predictions in some toy models 
[161, 162] and a general parameterisation can be provided in 
the framework of effective field theory [159]. For more details 
see the introduction by Ellis and Mavromatos [163] in Part A 
of this Special Issue. Typically, two scenarios are envisaged 
according to whether the linear term or the quadratic term is 
dominant, parametrised by the scale parameters %\ (linear case) 
and £2 (quadratic case) respectively. The point is that while 
QG effects would be prominent only at the Planck scales, there 
would be residual Lorentz Invariance Violation (LIV) effects 
at lower energies (GeV-TeV) in the form of anomalous photon 
velocity dispersion. 

Amelino-Camelia et al. [160] also noted that over a cosmo- 
logical distance L, the magnitude of time-delay Af induced by 
LIV between two photons with an energy difference AE is de- 
tectable: 



where a = 1 or 2 according to whether the linear or quadratic 
terms dominates in Eq. (3.1). The energy scale of QG is com- 
monly expected to lie somewhere within a factor g a of Ep. 
The best limit on the linear term has recently been placed by 
Fermi-LAT observations of GeV photons from GRB 090510 
(z = 0.903) which require M x = En/fr > 1.5 x 10 19 GeV 
[164]. The most constraining limit on the quadratic term M2 = 
Epi/fr > 6.4 x 10 10 GeV come from observations of an excep- 
tional flare of the active galactic nucleus (AGN) PKS 2155-304 
with the H.E.S.S. telescope [165]. 

It is important to keep in mind that although the QG induced 
time-delay is proportional to energy (as opposed to conven- 
tional dispersion effects which vary as inverse power of energy) 
similar time-delay effects may be intrinsic to the source [166]. 
Therefore, in order to distinguish between source and propaga- 
tion time-delays, different types of sources should be consid- 
ered with different physical properties and situated at different 
cosmological distances. For such studies, AGN and Gamma- 
Ray Bursts (GRBs) are the best candidates to test Eq. (3.2). 
AGN cover the higher energies (up to few TeV) and lower red- 
shift regime (probably up to z ~ 0.8 - 1) and GRBs the lower 
energies (probably few tens of GeV) but higher redshifts. Other 
promising candidates could be pulsars which until now have 
yielded constraints one order of magnitude weaker than the 
ones derived from AGN [167]. 

The consequences of improved sensitivity and larger energy 
coverage of CTA on time-delays recovery 

Using the Maximum Likelihood Estimation method (MLE) 
of Martinez and Errando [168], we investigate the effects that 
the improved CTA performance, in terms of increased statistics 
and broader energy lever arm, have on the time-delay recovery. 
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Five hundred realisations of Gaussian-shaped "pulsed" light 
curves were generated for several values of time-delays be- 
tween -60 sTeV" 1 and +60 sTeV" 1 in steps of 10 sTeV" 1 . 
This allowed an estimate of the value of the error 5t r on the mea- 

— 1/2 

sured time-delay Af r . The error decreases as N y , where N Y 
is the number of photons included in the likelihood fit, and sat- 
urates at a value of about 3 s TeV , about a factor of 3 less than 
the current generation of IACTs, due to the increased statistics 
of CTA. The effect of the increase in the energy lever-arm, pro- 
vided by the wide coverage of CTA from few tens of GeV to 
several tens of TeV, has also been addressed for the different ar- 
ray configurations taking into account the absorption of Extra- 
galactic Background Light (EBL) using the model of Kneiske 
et al. [169] 7 , and a spectral break at around 100 GeV. The pho- 
tons are thus separated into two populations at "high" and "low" 
energies with average values marked as Ehe and E^e respec- 
tively. Fig. 3.1 shows the variation of the energy lever-arm for 
different CTA configurations in energy-squared (quadratic case) 
AE 2 = E 2 HE - E 2 LE , after a convolution with the effective area 
and for a given choice of the energy value separating high and 
low energy bands Elim = 400 GeV. This value is almost stable 
regardless of the array or the spectral index [170]. The ranking 
shows that arrays /, C, J and H for the southern site and NB for 
the northern site are favourable for LIV studies. 

The intrinsic variability of the photon emission by astrophys- 
ical sources such as GRBs and AGN is the main systematic un- 
certainty in LIV searches. Until now, the detected variability of 
the AGN was limited to about 100 s, partially due to the limited 
statistics of the data. The possibility of improved separation 
of the initially unresolved double peak structures was investi- 
gated with light curve simulations and time-delay reconstruc- 
tion using again the MLE method. Fig. 3.2 shows the minimal 
peak separation which would allow distinguishing between two 
Gaussian spikes of the same standard deviation <x c for differ- 
ent photon statistics: a H.E.S.S.-like measurement, a CTA-like 
measurement with an improved photon collection by a factor 
100, and a more optimistic scenario with a factor 1000 more 
photons. 

Sensitivity to TeV photons for selected benchmark AGN 

Following the suggestion of Amelino-Camelia et al. [160] 
we define a sensitivity factor r\ = \At^\/T, where Afj is the 
magnitude of the time-delay introduced into the flare and T 
is duration of the burst/flare feature that is being examined. In 
Refs. [171, 172] it is shown that r\ > 0.3 is required to determine 
whether the observed flare time sequence has been skewed in 
comparison to its original form. To improve upon current lim- 
its will require observations of photons at energies larger than 
10 TeV (for a given redshift) or observations of similar flares 
from much more distant AGN. In the following, we calculate 
the integral numbers of photons from representative AGN to 
test LIV signatures in AGN flares. 

We have taken three VHE AGN representative of sev- 
eral known situations: an AGN-flare with high brightness 



7 This is a conservative choice — by using more transparent EBL models 
like that of Dommguez et al. [150], the results would be more promising. 
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Figure 3.1: Parameter AE 2 as a function of the spectral index for all arrays 
considered in CTA Monte Carlo simulations. The configurations have different 
layouts and number of telescopes. See details in the text. 
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Figure 3.2: Minimal distinguishable separation of two Gaussian peaks as a 
function of their width, in case of a H.E.S.S.-like measurement (open circles), a 
CTA-like measurement (crosses) and an ideal situtation where a ten times more 
photons than in the CTA measurement could be collected. 

(Mrk421), an AGN-flare that shows short variability timescales 
(PKS 2155-304) and the AGN with the largest known red- 
shift (3C 279) observed at VHE. The spectra in their highest 
recorded flux state have been taken from current IACT ob- 
servations, extrapolated to higher energies, convolved with the 
performance curves of the various CTA array layouts and inte- 
grated assuming a flare duration lasting for the appropriate time 
such that t] = 0.3. The results are shown in Fig. 3.3. Since this 
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falls into the category of unbinned methodologies, the precision 
in the time resolution is the same as the time precision of each 
array (i.e. better than 1 The uncertainty in the time-delay 
for a single photon is the time precision modulo of the energy 
resolution (being less than 10% as specified) and the distance 
(negligible), and saturates at 10 sTeV~' Gpc -1 . This translates 
into a Planck scale effect of (9(10) s (i.e. as good as a binned 
method would have for more than 100 photons). 

Mrk421 is known to show spectral hardening with increasing 
flux, and can have very hard spectra indeed on short timescales, 
as evidenced by Flare C in Ref. [173] — see discussion by Gai- 
dos et al. [174]. For the redshift z = 0.03 of Mrk421, Planck 
scale effects could be expected to induce a delay of ~ 1 s TeV~' ; 
for E > lOTeV photons this means we would need to be able 
to time resolve flare features at an unprecedented ~ 30 s du- 
ration. Whilst features this fast have yet to be identified, this 
could be because they are below the sensitivity of current in- 
struments and the top panel of Fig. 3.3 demonstrates that, if 
present, such features can indeed be probed. For PKS 2155- 
304, the redshift z = 0.117 implies that CTA would need fea- 
tures on the timescale of 120 s to test Planck scale effects, which 
is still a factor of a few faster than the 240 - 610 s rising and 
falling timescales of the ~1 Crab flares observed to date [175] 
but, as shown in the middle panel of Fig. 3.3, we would easily 
have sufficient photons to resolve such features. For 3C 279 
(z = 0.536, [176]), even though the flare timescale of 610 s re- 
quired for such a distant AGN are well within the variability 
timescales we currently observe for blazars, the photon flux we 
expect (~ 10% Crab) from such a distant source is expected to 
be too low to resolve such features at the highest energies, be- 
cause of the attenuation of the photon flux through interactions 
with the EBL. 

Concerning the different CTA arrays, Fig. 3.3 shows that the 
best performing arrays at high energies are C, D, H, I, J, K. 
While it may be of more interest to find out if they would detect 
sufficient photons on which to perform tests for time-delay, we 
note that there are a number of unbinned methods that can cope 
with sparse datasets [see, e.g., 171, 165, 177] so that 10 photons 
of E>=10 TeV are required in order to be able to begin to test 
for LIV. 

Time-delay recovery with realistic source lightcurve and spec- 
tra for the linear case 

For this study, we produce 2500 pairs of lightcurves (with 
and without time-delays), following the method of Timmer and 
Koenig [178] for each of the total 13 CTA arrays (11 South- 
ern array A. . . K and 2 Northern arrays NA, NB), as shown in 
Fig. 3.4. We scan the space of possibilities by selecting random 
values in the 5-dimensional space characterised by the follow- 
ing parameters: 

- time-delays in the range £i =0.1-5 (linear case) 

- AGN redshift linearly in the range 0.03 - 0.6 

- energy spectrum power law slope between 20 GeV and 
20 TeV in the range 1 - 2.5, and with a spectral cutoff at 
120 GeV 
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Figure 3.3: Integral number of events above a given energy, expected for the 
various array configurations for simulated flares, for N(> 10 TeV) = 10 pho- 
tons. Each panel is accumulated for the appropriate flare timescale required to 
be able to determine if Planck scale quantum gravity induced LIV is present. 
The top panel is for a 10 Crab Mrk421 flare [173] if it lasted for 30s duration; 
the middle panel is PKS 2155-304 similarly at its high level [175], but 120 s 
duration; the bottom panel is for a 610s 3C 279 flare at its highest recorded flux 
level [176]. 

- Flux level in the range 10~ n - 10~ 12 ph cirT 2 s 

- different observational periods: a) single day observations 
consisting of 3 pointings of 30 and 15 min, b) weekly ob- 
servations consisting of 3 nightly pointings of 30 and 15 
min, and c) monthly observations consisting of 2 nightly 
pointings of 30 and 15 min. 

The photons thus generated are then distributed as a function 
of time, based on the variability type that we have initially as- 
sumed (e.g. red-noise). These light curve pairs incorporate de- 
lay effects accumulated over a given distance depending on the 
AE of each pair. The light curve pairs are subsequently con- 
volved with the CTA arrays performance, using the effective 
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area, the background count rate and differential sensitivity. Fi- 
nally, for each CTA array, we recover the observed time-delays 
using the cross-power spectral analysis method of Nowak et al. 
[179]. For each pair of light curves and for each array we con- 
sider the quality factor between the simulated time-delay Afd 
and the recovered time-delay Af r defined as: 

q = \At d - At r \/6t r , (3.3) 

where 6t r is the width of the Gaussian distribution of the recov- 
ered time-delays coming from the simulations. 
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Figure 3.4: Time-delay recovery fraction with quality factor q < 1 for each 
CTA configuration. See text for details. 

When considering all the 13 arrays together, we find on aver- 
age that 6% of the time-delays are recovered with q < 1. This 
is a very strict limit. If one relaxes this limit to q < 2 and q < 3, 
respectively 77% and 99% of the time-delays in our sample of 
events are recovered, thus making the prospects of detecting (or 
constraining) LIV signatures with CTA rather optimistic. 

To understand which arrays have the best prospects, the 
time-delay recovery results for each array individually is shown 
in Fig. 3.4, adopting the limit of q < 1. From the plot we see 
that the best CTA configurations for time-delay recovery are 
C, D, H, I and NB respectively for the Southern and Northern 
hemispheres. 

4. Other physics searches with CTA 

In this section, we highlight topics of fundamental physics 
searches that were discussed in recent years and whose scenar- 
ios could be studied or hopefully constrained with CTA. 

There are several caveats if one wants to present such various 
and complex topics "in a nutshell". First of all, the list of topics 
is not exhaustive; only a subset of topics is reported here. Sec- 
ond, some of the studies presented in this section may not be 



up-to-date by the time this article is published: theories in this 
area evolve and are updated exceedingly rapidly. In addition, 
most of these studies were formulated only within the context 
of the current generation of IACTs, and not for CTA. When- 
ever possible, considerations about the prospects for CTA will 
be addressed. Third, the discussion will mostly be of only a 
qualitative nature. The goal in this contribution is to provide an 
introductory discussion of the area, with the aim of encouraging 
others to explore in more detail these and other interesting new 
physics possibilities. Let us add that pursuing exotic physics 
with IACTs (and hence CTA) should be done because a) it is 
possible (this may seem a naive argument, but, given the terra 
incognita offered by a new observatory such as CTA, it is a 
strong one; b) VHE gamma rays have been identified as likely 
drivers of truly fundamental discovery. VHE gamma rays are 
a tool to explore new physics and new astrophysical scenarios, 
the nature of which may contain yet unknown, and unexpected, 
features. The potential for revolutionary discovery is enormous. 

4.1. High energy tau-neutrino searches 

Although optimized to detect electromagnetic air showers 
produced by cosmic gamma-rays, IACTs are also sensitive 
to hadronic showers. Inspired by calculations made by the 
AUGER collaboration and D. Fargion [see, e.g., 180, 181, 182, 
183, 184, 185], the possible response of IACTs to showers ini- 
tiated by very high energy r-particles originating from a v T col- 
lision with the sea or underneath rock is described. 

It is well known that neutrinos of energies above the TeV 
energy range can form part of the cosmic rays hitting the 
Earth. The origin of such neutrinos could be from point-like 
sources like galactic microquasars [186, 187] or extragalactic 
blazars [188, 189] or gamma-ray bursts [190]. There are also 
diffuse fluxes of high energy neutrinos predicted to come from 
unresolved sources, including interactions of EHE cosmic rays 
during their propagation [191]. Finally, one could think of a 
more exotic origin of high energy neutrinos like those coming 
from DM particle annihilation, topological defects or cosmic 
strings [192, 193, 194]. Neutrinos are produced in astrophysi- 
cal sources or during the transport, mainly after pion and sub- 
sequent muon decays: 

7T + -> // + v„ yt -> v„ e + Ve (4.1) 

n~ -> iT v M iT -> Vfj e~ v e (4.2) 

such that the typical neutrino family mixing at the source is 
(v e ,Vfj,v T ) = (1:2:0). Tau-neutrinos are found either at the 
source, if charmed mesons are formed instead of pions, or are 
created during the propagation, after flavor mixing, such that at 
Earth, the neutrino family mixing could be (v e , v^, y T )Enrth = 
(1:1:1) [195]. 

The v T channel has several advantages with respect to the 
electron or muon channel. First, the majority of the possible 
t decay modes lead to an (observable) air shower or a combi- 
nation of showers. Only 17.4% of the decays lead to a muon 
and neutrinos, considered to be unobservable for the effective 
areas of interest here. Moreover, the boosted t lifetime ranges 
from some 50 m at 1 PeV to several tens of kilometers at EeV 
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energies, almost unaffected by energy losses in matter and thus 
surpassing the muon range by a factor of ~ 20. Finally, the 
originating r decays, instead of being absorbed by matter, and 
thus gives origin to another v T of lower energy which in turn can 
produce a r. At the highest energies, the Earth becomes com- 
pletely opaque to all types of neutrinos giving rise to a pile-up 
of v T - T. 

To be able to observe atmospheric showers from v T , the tele- 
scopes should be pointed at the direction where the r escapes 
from the Earth crust after having crossed an optimized distance 
inside the Earth. Of course, this distance is strongly dependent 
on the telescope location, and no general conclusions can be 
drawn before the CTA site will be defined. In the past, two 
directions were proposed 8 : the observation slightly below the 
horizon downhill, e.g. if the telescope is located at a moun- 
tain, and the observation through a possible mountain chain in 
the vicinity (see Fig. 4.1). Both observations are at extremely 
low elevation angles, i.e. the telescopes pointing horizontally or 
even below the horizon, a fact that guarantees that the hadronic 
background diminishes until almost vanishing at the horizon. 
Only a small light contamination from continuous scattered star 
light and from scattered Cherenkov light by air showers will 
then be observed. 




Figure 4.1: Sketch of v T searches with CTA. CTA telescopes can be pointed 
horizontally to a far mountain (if present) and observed the emerging T-induced 
atmospheric shower. 

In Ref. [196], the effective area for v T observation with the 
MAGIC telescope was calculated analytically. The results were 
the following: the maximum sensitivity would be in the range 
100 TeV-1 EeV. For the observation downward towards the Sea, 
the sensitivity for diffuse neutrinos is very low because of the 
limited FOV (3 events/year/sr) and CTA cannot be competi- 
tive with other experiments like Icecube [197], Baikal [198], 
Auger [184], Antares [199] or KM3NeT. On the other hand, 
if flaring or disrupting point sources are observed, like is the 
case for GRBs, one can even expect an observable number of 
events from one GRB at reasonable distances, if the event oc- 
curs just inside a small observability band of about 1 degree 
width in zenith and an azimuth angle which allows to point the 
telescopes downhill. 

For CTA, the situation could be different: taking an extension 
of the FOV of several times that of MAGIC in extended obser- 
vation mode, the higher effective area and lower energy thresh- 
old, meaning higher fluxes, one naive rescaling of the MAGIC 

8 The study was made for the case of MAGIC [196] 



calculations leads to relatively optimistic results, depending 
very much on the local geography. For point-like sources, the 
situation would not change so much w.r.t. the MAGIC case, un- 
less the CTA telescopes are located close to a shielding moun- 
tain chain. The required observation times are still large, but 
one may argue that these observations can be performed each 
time when high clouds preclude the observation of gamma-ray 
sources. 

4. 2. Ultrarelativistic Magnetic Monopoles 

The existence of magnetic monopoles is predicted by a 
wide class of extensions of the standard model of particle 
physics [200]. Considerable experimental effort has been 
undertaken during the last eight decades to detect magnetic 
monopoles. No confirmed success in detection has been re- 
ported at the present time. Current flux limits on cosmogenic 
magnetic monopoles reach values of O(10~ 15 cm _2 s _1 sr _1 ) to 
<9(10~ 17 cm _2 s _1 sr _1 ) depending on the monopole velocity. As 
outlined below, the CTA observatory is sensitive to a magnetic 
monopole flux. 

According to Tompkins [201] magnetic monopoles moving 
in air faster than the speed of light in air are emitting * 4700 
times more Cherenkov photons than an electric charge under 
the same circumstances. Being fast enough (Lorentz factor 
y > 10 3 ) and heavy enough (mass Mc 2 > 1 TeV) magnetic 
monopoles that possibly propagate through the earth atmo- 
sphere are neither significantly deflected by the Earth's mag- 
netic field nor loose a significant amount of energy through 
ionization [202]. Assuming the last two constraints to be ful- 
filled a magnetic monopole moving through the Earth's atmo- 
sphere propagates on a straight line, thereby emitting a large 
amount of Cherenkov photons. This process of a uniform emis- 
sion of intensive Cherenkov light differs from the Cherenkov 
light emitted by secondary particles in a shower initiated from 
a high energy cosmic or gamma-ray. As shown by Spengler 
and Schwanke [202], the number of triggered pixels in a tele- 
scope array is typically smaller and the intensity of the triggered 
pixels is typically higher for magnetic monopoles compared to 
events originating from cosmic or gamma-rays. Cuts in a pa- 
rameter space spanned by the number of triggered pixels in the 
CTA array and the number of pixel with high intensity allow for 
an excellent discrimination between magnetic monopole events 
and background from cosmic or gamma-rays. The effective de- 
tection area of H.E.S.S. [203] for magnetic monopoles has been 
studied in detail [202]. Extrapolating the results of this study for 
CTA with its one order of magnitude increased design collec- 
tion area, leads to a typical CTA magnetic monopole effective 
area of 4500 m 2 sr. In Fig. 4.2, we show that assuming around 
3000 hours of CTA data from different observations accumu- 
lated in about 4 years of array operation, the sensitivity of CTA 
to magnetic monopoles with velocities close to the speed of 
light can reach the Parker limit [200] of O(10 15 cin -2 s _1 sr _1 ). 
Despite being still two orders of magnitude worse than current 
monopole flux limits from neutrino experiments [204] this sen- 
sitivity will allow a technically independent and new test for the 
existence of magnetic monopoles. 
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Figure 4.2: Upper limits on the magnetic monopole flux as function of the 
monopole velocity inferred from different experiments [200, 204] together with 
the estimated sensitivity of CTA. The sensitivity range shown for CTA indicates 
the dependence of the sensitivity on CTA design parameters and analysis de- 
tails. 

4.3. Gravitational waves 

The period of operation of CTA should hopefully see the de- 
tection of the first gravitational wave (GW) by ground-based in- 
terferometers, now in the "advanced sensitivity" design phase. 
The 3km-scale Michelson interferometers Enhanced LIGO 
and Advanced Virgo [205] are increasing their sensitivity and 
extending the horizon distance of detectable sources up to hun- 
dreds of Mpc, depending on the frequency. Two additional 
smaller interferometric detectors are part of the network of 
GW observatories, the Japanese TAMA (300 m arms) and the 
German-British 600 m interferometer GEO600. The LIGO and 
Virgo observatories should start full operation in the advanced 
version in 2014/2015 and may operate together to fully recon- 
struct the arrival direction of a signal. They may localize strong 
GW bursts with an angular uncertainty down to one degree, 
while weaker signal have larger uncertainties, up to tens of de- 
grees [206]. At the stage of conceptual design, the Einstein 
Telescope aims at increasing the arm length to 10 km, with three 
arms in a triangular pattern, implementing consolidated tech- 
nology. It is foreseen to be located underground to reduce the 
seismic motion thus allowing a better sensitivity up to a factor 
10 [207]. 

The most promising astrophysical mechanisms able to pro- 
duce observable GWs are in-spiral and coalescence of binary 
compact objects (neutron stars and black holes), occurring for 
example during the merger of compact binary systems, as pro- 
genitors of supernova or neutron-star collapse, and associated 
with pulsar glitches. Signals from these systems may last from 
milliseconds to a few tens of seconds, but their expected rate 
and strength are uncertain [for a review, see, e.g., 208]. More- 
over, unexpected or unknown classes of sources and transient 
phenomena may be responsible for GW emission and may ac- 
tually provide the first detection. Therefore combined GW and 
electromagnetic observations would be critical in establishing 
the nature of the first GW detection. While electromagnetic 
counterparts cannot be guaranteed for all GW transients, they 
may be expected for some of them [209, 210] from radio waves 



to gamma-rays, such as in gamma ray bursts [211], ultra high- 
luminous X-ray transients and soft gamma repeater [212]. Elec- 
tromagnetic identification of a GW would confirm the GW de- 
tection and improve the reconstruction and modeling of the 
physical mechanism producing the event. Moreover, significant 
flaring episodes identified in the electromagnetic band could 
serve as an external trigger for GW signal identification, and 
could even be used to reconstruct independently the source po- 
sition and time, thus allowing the signal-to-noise ratio required 
for a confident detection to be lowered. The feasibility of this 
approach has been corroborated through dedicated simulations 
by the LIGO and Virgo collaborations [213]. 

CTA has the capability to pursue such a program of immedi- 
ate follow-up of target of opportunity alerts from GW observa- 
tories, and to interact with GW collaborations to pursue offline 
analysis on promising candidates. The capability of CTA to 
observe in pointed mode with small FOV or in extended mode 
covering many square degrees of sky, is unique to follow strong 
and weak GW alerts. The observation mode should resemble 
the GRB procedure, which allow a fast repositioning on the or- 
der of tens of seconds. 

In the era in which ground-based gravitational wave detectors 
are approaching their advanced configuration, the simultaneous 
operation of facilities like CTA and Virgo/LIGO may open, in 
the forthcoming years, a unique opportunity for this kind of 
multi-messenger search. 

5. Summary and Conclusion 

In this study we have investigated the prospects for detection 
and characterization of several flavors of physics beyond the 
standard model with CTA. 

Particle Dark Matter searches 

We have investigated dark matter (DM) searches with CTA 
for different observational strategies: from dwarf satellite 
galaxies (dSphs) in Section 1.1, from clusters of galaxies in 
Section 1 .2 and from the vicinity of the Galactic Centre in Sec- 
tion 1.3. In Section 1.4, we discussed spatial signatures of DM 
in the diffuse extragalactic gamma-ray background. 

Concerning searches in dSphs of the Milky Way, we have in- 
vestigated the prospects for detection of well-known "classical" 
dSph like Ursa Minor and Sculptor, and one of the most promis- 
ing "ultra-faint" dSph, Segue 1 (Table 1.1). We have first shown 
that the predictions for core or cusp DM density profiles are 
quite similar for the baseline CTA angular resolution (Fig. 1.1). 
We have then simulated a 100 h observation for several CTA 
arrays, and found that for Segue 1, we can exclude velocity- 
averaged cross-sections (cr aml v) above 10~ 23 - 10~ 24 cm 3 s _1 de- 
pending on different annihilation channels (Fig. 1.2). We also 
presented the same results in terms of the minimum astrophys- 
ical factor for dSphs to be detected (Fig. 1.3), showing that as- 
trophysical factors of at least 10 21 GeV 2 cirT 5 are needed. We 
finally showed the minimum intrinsic boost factor to achieve 
detection (Fig. 1 .4), which for Segue 1 is about 25 for a hard an- 
nihilation spectrum. The best candidate arrays for dSph study 
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are array B and E. Nevertheless, the robustness of our results 
is hindered by the yet not precise determination of the astro- 
physical factor in some cases. Forthcoming detailed astronom- 
ical measurements will provide clues for deep exposure obser- 
vations on the most promising dSphs, with, e.g., the planned 
SkyMapper Southern Sky Survey [214], which will very likely 
provide the community with a new dSph population, comple- 
menting the Northern hemisphere population discovered by the 
SDSS. Also, the uncertainties on dark matter density will be 
significantly reduced by new measurements of individual stel- 
lar velocities available after the launch of the GAIA mission 9 . 
Stacking-methods of Fermi-LAT dSphs data were proven valid 
to make constraints more stringent [49, 215, 216]. The applica- 
tion of these methods for CTA is currently under study. 

The search for DM signatures in galaxy clusters, investigated 
in Section 1.2 was performed for two representative clusters, 
Perseus and Fornax. The former one is thought to have the 
highest CR-induced photon yield, and the latter is thought to 
have the strongest DM-induced signatures. Compared to dSphs, 
the gamma-ray signatures of galaxy clusters have several con- 
tributions: in the first place, the DM signal is expected from 
an extended region that can be larger than a few degrees, and 
secondly, gamma-rays induced by interactions of accelerated 
cosmic rays with the ambient fields and/or by individual clus- 
ter galaxies are an irreducible background to the DM signal, 
as recently shown in Refs. [103, 102]. We have simulated the 
prospects of detection in 100 h of observation by using MC 
simulations of extended sources. Regarding DM signatures, we 
have used the model of Pinzke et al. [30] for the Fornax cluster, 
and showed that in 100 h we could put contraints on the order of 
(cr ann v) < 10~ 25 cm 3 s~' (Fig. 1.6), which are competitive with 
respect to those obtained with dSphs. The results are promis- 
ing: if the intrinsic boost factor from subhalos is larger than 
that predicted by the model we used, or mechanisms of Som- 
merfeld enhancement are at work, there is also the possibility 
to have a detection in 100 - 200 h with array B or E. We have 
also considered the prospects of detection of CR-induced sig- 
nal in hadronic acceleration scenarios in Fig. 1.5. We have seen 
that the CR-induced emission from the Perseus cluster could be 
detected in about 100 h. Finally, we discussed the more real- 
istic case when DM- and CR-induced gamma-rays are treated 
together. We discuss that the difference in both the spatial and 
spectral features of the two emissions can be used as a method 
for discrimination, while more quantitative results need dedi- 
cated MC which were not available when writing this contribu- 
tion. We underline that the extension of the expected DM emit- 
ting region in galaxy clusters represents a problem for current 
Cherenkov Telescopes since their FOV is limited to 3 - 5 de- 
grees and their sensitivity rapidly decreases moving away from 
the centre of the camera. CTA will overcome this limitation, 
having a FOV of up to 10 deg and an almost fiat sensitivity up 
to several degrees from the centre of the camera. For galaxy 
cluster searches, CTA will hence mark the difference compared 
to the current generation of IACTs. 
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More promising are DM searches of annihilation signatures 
in the Galactic halo, where the DM density is expected to be 
known with much higher precision than in the Galactic Cen- 
tre itself or in (ultra-faint) dSphs or galaxy clusters. This was 
studied in Section 1.3. By adopting dedicated observational 
strategies of the region close to the Galactic Centre, as shown 
in Fig. 1.8, it was shown that CTA has the potential to reach the 
thermal annihilation cross-section expected from WIMP DM 
of 10~ 26 cm 3 s~' and lower (Fig. 1.9) in 100 h observation of the 
vicinities the Galactic Centre using the "Ring" method. Models 
with a large photon yield from DM annihilation will be con- 
strained for even smaller cross-sections. It is also expected that 
the limits presented here can be improved by factor of a few 
when the stereoscopic analysis of CTA events has been under- 
stood so well that a further suppression of the background be- 
comes feasible. This would be the first time that ground-based 
Cherenkov telescopes could reach this sensitivity level. 

Besides observations of individual dedicated objects, the 
capabilities of CTA for searching DM signals in the diffuse 
background of gamma-ray radiation were discussed in Sec- 
tion 1.4. We discussed the reconstruction performance for 
different anisotropy power spectra and residual background 
level. Considering a current model for the anisotropy power 
spectra, we showed that CTA may be able to distinguish a DM- 
induced diffuse gamma-ray component from the astrophysical 
background. 
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Figure 5.1: Comparison of exclusion curves of Fermi-LAT in 24 months [52] 
and expected for 10 years (rescaled with the square root of time). The exclu- 
sion curves for the various targets studied in this contribution are also reported 
for the bb annihilation channel: for the dwarf satellite galaxy Segue 1 (green 
curve, see Sec. 1.1), for the Fornax galaxy cluster in case only DM-induced 
gamma-rays are considered (blue line, see Sec. 1.2) and for the ring-method of 
observation of the Galactic Centre vicinities (red line, see Sec. 1.3). 

In Fig. 5.1, we summarize the constraints that we expect 
with CTA for a WIMP annihilating purely into bb in 100 h 
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observation, with the different targets discussed above. As al- 
ready anticipated, the best results are expected for the observa- 
tion of the vicinity of the Galactic Centre, where we expect to 
reach the thermal annihilation cross-section for WIMP DM of 
10~ 26 cm 3 s~'. Unlike present IACTs, whose sensitivity super- 
sedes that of the Fermi-LAT at masses around a TeV, CTA will 
constitute the most sensitive instrument above masses of about 
100 GeV. It should be noted that these estimates are conser- 
vative: the most important improvement can be expected from 
the possible redefinition and final optimization of the array lay- 
out 10 In addition, the presented sensitivities were calculated us- 
ing generic analysis which was not optimized specifically for 
the DM searches and thus our results could be considered con- 
servative in this sense. 

Obviously, a firm identification of DM requires a very 
good spectral discrimination with respect to any possible 
astrophysical background. Spectral shapes and (even more 
so) absolute normalization of these backgrounds are often 
poorly determined and the DM signal most likely is small 
in comparison. As the extent to which these factors affect 
detection claims is highly model and target dependent, we refer 
to future detailed more focused assessments. 

It has been shown that the detection of gamma-rays pro- 
vides complementary information to other experimental probes 
of particle DM, especially that of direct detection, because 
CTA could be able to access a fraction of the parameter space 
not accessible otherwise [217, 218]. With respect to particle 
searches at the LHC, the comparison is not straightforward, as 
LHC results are usually strongly related to specific models, and 
general conclusions are somewhat model dependent, as shown 
by recent publications from the ATLAS and CMS collabora- 
tions [219, 220, 221]. Generically, the discovery of a candidate 
for particle DM will be limited by the available centre-of-mass 
energy. Other scenarios exist, in the context of specific super- 
symmetric models for DM, that exhibit parts of the model space 
not accessible by the LHC [222]. In any case, LHC discovery 
of dark matter, would prompt the need for proof that the particle 
is actually consistent with the astrophysical DM, and close col- 
laboration with LHC physicists is currently under organization 
to facilitate the optimal use of accelerator results within CTA. 
A concrete scenario has been analyzed by Bertone et al. [223] 
in the case of a SUSY model in the so-called co-annihilation 
region. Simulated LHC data were used to derive constraints on 
the particle physics nature of the DM, with the result that the 
LHC alone is not able to reconstruct the neutralino composi- 
tion. The situation improves if the information from a detection 
of gamma-rays after the observation of the Draco dSph by CTA 
is added to the game: in this case the internal degeneracies of 
the SUSY parameter space are broken and including CTA al- 
lows us to fully interpret the particle detected at the LHC as the 
cosmological DM. In the other case where the LHC will not de- 



10 For example, it is currently under disussion, the possibility to add 36 
medium-size telescopes of Schwarzchild-Couder design to the arrays consid- 
ered here. Preliminary simulations predict that improvement in the overall sen- 
sitivity by at least a factor of 2 compared to that studied here may be expected. 



tect any physics beyond the Standard Model, predictions were 
made in the context of the CMSSM [224] indicating that the 
mass of the neutralino will be bound to be larger than approxi- 
mately 250 GeV (400 GeV) if any new physics will be detected 
by the LHC for an energy of the centre-of-mass yfs = 14 TeV 
and a luminosity of 1 fbr 1 (100 fbr 1 ). In this scenario, CTA 
could be the only instrument to be able to detect and identify a 
WIMP candidate with masses beyond some hundreds GeV. 

Axion-like particle searches 

In Section 2, the prospect of searches for axion-like parti- 
cle (ALP) signatures with CTA were studied. We saw that the 
theoretical photon/ ALP mixing has important implications for 
astronomical observations, in such that the mixing could dis- 
tort the spectra of gamma-ray sources, such as Active Galactic 
Nuclei (AGN) (or galactic sources), in the TeV range. This dis- 
tortion adds to that caused by the absorption of the gamma-ray 
photons with UV and IR photons of the Extragalactic Back- 
ground Light (see Fig. 2.1). The photon flux recently measured 
by some experiments, in particular at TeV energies, already 
exceeds that predicted by conventional models which attempt 
to explain spectra in terms of observed source spectra and/ or 
EBL density [144, 145, 146, 147, 148], one should not expect a 
photon flux as high as recently measured by some experiments, 
in particular at TeV energies. The hard spectrum deduced for 
some AGN is difficult to explain with conventional physics as 
well. While it is still possible to solve these puzzles without 
exotic physics, photon/ ALP conversions may naturally allevi- 
ate both problems. In order to quantitatively study the effect of 
photon-axion conversion over cosmological distances, the to- 
tal photon flux from a simulated flare of a far-distant source 
was considered. The source was simulated based on the flat 
spectrum radio quasar 4C +21.35 (PKS 1222+21, z = 0.432, 
based on the observation performed by MAGIC [149]), assum- 
ing an intrinsic unbroken power-law spectrum, and trying to 
understand the observability under different absolute flux nor- 
malization and flare duration (Figs. 2.2 and 2.3). The range of 
characteristic scale energy (critical energy, E crit ) described in 
Eq. (2.3), and thus the ALP mass that can be probed with CTA 
for the different ALP scenarios, is unknown and may be tested 
with CTA. In general, the distortion of the spectra due to ALP 
depends on the particular case, but as a general trend it will 
become larger as we increase the observation time and/or the 
flux of the source. As an example, we found that a 0.5 h dura- 
tion flare like the one reported by MAGIC would not be enough 
for CTA to detect a significant effect in any of the tested ALP 
scenarios (Fig. 2.4). However, a flare of similar intensity, but 
lasting 5 h would already be enough to see the boost due to 
ALPs for those scenarios with E en( < 500 GeV. For a hypothet- 
ical flare with an intensity 5 times larger, lasting 5 hours, the 
accessible range of E en( would extend up to 1.3 TeV (Fig. 2.4). 
Hopefully, not only PKS 1222+21 but also many other similar 
objects will be followed-up by CTA in the near future, making 
the field of ALP searches very promising. 

We must emphasize that a boost in the flux is only possible in 
the energy range where the EBL is already at work. Thus, even 
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for the most distant sources detected to date by IACTs, the en- 
ergy range below ~ 100 GeV would not probably be of much 
help. On the other hand, even when E cr i, lies within the energy 
range covered by IACTs, the drop/jump might not be accessi- 
ble to these instruments. This would be the case, for instance, if 
E cr n is at the highest energies, from several to tens of TeV: the 
attenuation due to the EBL for a distant source would be huge, 
and the resulting flux, even after accounting for the ALP boost, 
too low to be detected by current IACTs or by CTA under any 
of the possible array configurations. Taking all these consider- 
ations into account, the most suitable energy for ALP searches 
with CTA seems to be an intermediate one in which the EBL is 
already present but still introduces only a moderate absorption, 
i.e. from a hundred GeV to a few TeV. As a result, we do not 
expect to obtain largely different results with candidate array 
configurations other than the one we used (array E), since they 
all perform very similarly in the intermediate energy range. 

Finally, although very challenging given the uncertainty in 
the value of the IGMF, we should stress that the lack of detec- 
tion of suspicious features in the spectra of distant gamma-ray 
sources might translate into useful constraints of the ALP pa- 
rameter space (coupling constant and ALP mass). A more de- 
tailed study is definitely needed in order to find out what should 
be the best strategy to achieve the strongest constraints. This 
study will be done elsewhere. 

Lorentz Invariance Violation 

In scenarios where Lorentz invariance is violated by quantum 
gravitational effects, the space-time fabric may be distorted so 
that the vacuum shows a non-unitary refractive index and thus 
the light speed would be wavelength dependent. Observation 
of gamma-ray flares from far distant objects like active galactic 
nuclei or gamma-ray bursts, may allow to detect the time-delay 
between photons of different energies not caused by intrinsic 
source mechanisms. In Section 3, we discussed the sensitivity 
of the different CTA array configurations on detecting time- 
delays induced by Lorentz Invariance Violations (LIV). While 
limits on LIV from the current generation of IACTs are weaker 
than those estimated from Fermi-LAT measurement in the so- 
called linear case, CTA is likely to invert this scenario. 

Using for the Maximum Likelihood Estimation method of 
Martinez and Errando [168], 500 Gaussian-shaped pulsed light- 
curves with time-delay from -60 s TeV -1 to 60 s TeV -1 were 
simulated and reconstructed with the different arrays. CTA will 
have improved statistics of photons and larger spectral lever- 
arm due to the enlarged energy range with respect to the cur- 
rent generation of IACTs. This will allow to better differentiate 
between the two Gaussian peaks as shown in Fig. 3. 1 for the dif- 
ferent arrays. This ability to differentiate peaks is also discussed 
in Fig. 3.2 for different width of the peaks. The best array con- 
figurations were discussed. As a result of these studies, a gain 
of about a factor 50 in the LIV scale for the "quadratic" model 
(see Eq. (3.1)) is expected, compared to current generation of 
telescopes, while the limits on the linear term will largely ex- 
ceed the Planck energy scale. 

We then used extrapolation to high-energies of real AGN 
spectra observed by the current generation of IACTs for three 



representative scenarios: a very bright AGN (Mrk 421), a fast- 
variable one (PKS 2155-304) and a high-redshift one (3C 279). 
High-energy photons above 10 TeV will guarantee the best sen- 
sitivity to observe LIV signatures, and CTA with its improved 
sensitivity at those high energies, will allow to collect sufficient 
photons, whereas photon statistics will always be the final lim- 
iting factor on tests for time-delay. 

Finally, pairs of realistic AGN lightcurves with and without 
time-delay were simulated, and folded with CTA performance. 
For each array, we calculated the fraction of photons in which 
the time-delay was successfully measured according to a qual- 
ity factor q (Eq. 3.3). In the most stringent case (q < 1), we 
report the photon fraction recovery of each individual array in 
Fig. 3.4. As a main result, more than 10% of the time-delays 
can be recovered with several possible CTA arrays. If we relax 
the quality factor, and thus the precision on the reconstructed 
time-delay, essentially all the time-delays are recovered. We 
showed that arrays C, D, H, I and NB, respectively for the 
Southern and Northern hemisphere, have the best chance for 
detection. 

Based on these genuinely different time-delay reconstruc- 
tion methods we ensure that our final results, with respect to 
CTA-array ranking, are free from any possible systematic ef- 
fects related to a given analysis method e.g. idealized source 
redshift-distribution, idealized source light-curves, and ideal- 
ized source time-scales. In all analyses methods the sub-arrays 
C, H, I and NB seem to be sufficiently good to perform de- 
tection of LIV effects by measuring differences in the arrival 
times of VHE photons. That means that these arrays are in gen- 
eral sufficiently good to perform temporal studies of light-curve 
signals and even detection of time-delays in AGN induced in- 
trinsically in the source. The latter is an interesting degeneracy, 
connected with the actual origin of the time-delays, that CTA 
will definitely be able to break through population studies not 
based on exceptional flaring states but on a routine basis. 



Other Searches 

Finally, in Section 4, we have qualitatively discussed 
the physics case of a selection other exotic physic searches 
which are in principle possible with CTA: the observation of 
atmospheric showers from T-particles emerging from the Earth 
crust, the observation of atmospheric showers from magnetic 
monopoles, and the possible follow-up of gravitational waves 
events. Despite the prospects being sometimes pessimistic, 
those subjects were shown to underline again the possibility of 
using an astronomical observatory such as CTA for fundamen- 
tal physics searches. 

As a final closing remark, we believe that CTA could offer 
one of the most powerful tools in the study of some of the most 
pressing questions in modern physics. In the next few years it 
may lead to a range of new observables, new methods and new 
theories. In preparation for these developments, it is essential 
that work such as that performed here is continued. 
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